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We develop a theory of parametric resonance in tunable superconducting cavities. The nonlin- 
earity introduced by the SQUID attached to the cavity, and damping due to connection of the 
cavity to a transmission line are taken into consideration. We study in detail the nonlinear classical 
dynamics of the cavity field below and above the parametric threshold for the degenerate para- 
metric resonance, featuring regimes of multistability and parametric radiation. We investigate the 
phase-sensitive amplification of external signals on resonance, as well as amplification of detuned 
signals, and relate the amplifier performance to that of linear parametric amplifiers. We also discuss 
applications of the device for dispersive qubit readout. Beyond the classical response of the cavity, 
we investigate small quantum fiuctuations around the amplified classical signals. We evaluate the 
noise power spectrum both for the internal field in the cavity and the output field. Other quantum 
statistical properties of the noise are addressed such as squeezing spectra, second order coherence, 
and two-mode entanglement. 

PACS numbers: 85.25.-j, 84.30.Le, 84.40. Dc, 42.50.Lc, 42.65.Yj 



I. INTRODUCTION 

Parametric resonance is a fundamental physical phe- 
nomenon that is encountered eventually in every area 
of science. In different disciplines, however, different 
facets of this rich phenomenon play a major role and 
are highlighted. Parametric instability and multistablc 
regimes in nonlinear dynamicJ^, noise d rive n transitions 
among stable states in statistical physicJ^El, wave mixing 
and frequency conversion in wave dynamic^ are topics 
of primary interest. In electrical and optical engineer- 
ing the low-noise properties of parametric amplifiers at- 
tract attentions, as well as non-classical statistical prop- 
erties of the electromagnetic field generated by paramet- 
ric device^^t^. 

In superconducting electronics, the idea of using 
Josephson junctions for quantum limited parametric am- 
plification is under attention and development since the 
During the last years the field revived by chal- 
lenges of quantum information technology. The circuit- 
QED design, ini tially proposed for qubit manipulation 
and measuremenlfi^E^, was employed for developing a va- 
riety of parametric device9l^"I21 

The circuit-QED approach is based on a combination 
of extended linear electromagnetic elements (transmis- 
sion lines and resonators) with Josephson junctions as 
nonlinear lumped elements. The design is flexible, allow- 
ing for diverse methods of parametric pumping, phase 
preserving and phase sensitive amplification schemes, dif- 
ferent numbers or inpu t and output ports, distributed 
Josephson nonlinearitie j^° I ^ H 

The most of developed amplifiers are engineered in 
such a way that the dominant pump tone is sent through 
the same port as the signal, and parametric resonance 
is achieved by mixing them in nonlinear Josephson ele- 
ments. A different method is available for tunable su- 
perconducting cavitie^22123j^ j^^ie device consists of a res- 
onator terminated with one (or more) dc-SQUID(s) that 



determines the reflection condition at the cavity edge and 
hence the cavity resonance spectrum. Parametric res- 
onance is achieved by rapid modulation of a magnetic 
flux through the SQUID with an appropriate frequency. 
A number of interesting parametric effects have been 
observed with such a device: phase sensitive ampliflca- 
tioni^, frequency conversioiP^, radiation and multistabil- 
ity regimes above the parametric thresholcP^, quantum 
entanglement of output photon^, generation of pho- 
tons out of v acuum noise^Zl - an analog of the dynamical 
Casimir effecPI^Il. 

In this paper we formulate a consistent theory of para- 
metric resonance in a tunable superconducting cavity. 
We aim at a unifled picture of the phenomenon below 
and above the parametric threshold. To this end we 
include into consideration the SQUID nonlinearity, and 
damping due to connection to a transmission line. The 
latter provides a stage for studying the parametric am- 
pliflcation. We develop a full nonlinear description of the 
cavity resonance dynamics and the ampliflcation effect in 
the classical limit, and study small quantum fluctuations 
of amplifled and radiative fields. For certainty we con- 
sider parametric excitation of the main cavity mode ujq 
by pumping with a frequency close to twice the cavity 
resonance, fl ~ 2ujq. 

The overall picture of nonlinear parametric resonance 
in the tunable cavity is rather rich and complicated. At 
very small pump strength the cavity intrinsic dynamics 
resembles the one of the Duffing oscillato/^ showing a bi- 
furcation of the cavity response and bistability. However, 
the scattering of an external incidental wave is qualita- 
tively different from the Duffing case: the scattering is 
inelastic, the reffected wave undergoes ampliflcation or 
deampliflcation depending on the phase shift between the 
input tone and the pump (phase sensitive ampliflcation) . 

With increasing pump strength, the amplification ef- 
fect increases, and at the same time the resonance nar- 
rows such that the bifurcation occurs at ever smaller in- 



put amplitudes. Eventually, while approaching the para- 
metric threshold, the cavity response becomes nonlinear 
at any small input amplitude. 

Further increase of the pump strength leads to an in- 
stability of the cavity zero-amplitude state and the for- 
mation of finite-amplitude states accompanied by sta- 
tionary parametric radiation at the half frequency of the 
pump. The radiative states are bistable in a certain win- 
dow of detuning of the pump frequency from the cavity 
resonance. Outside of this interval at red detuning the 
radiative states coexist with the stable zero-amplitude 
state (tristability) , and the latter one becomes dominant 
at far red detuning. Remarkably all these multistable 
regimes have been observed in experiment with a high 
quality tunable cavitjEH. 

The multistability regimes are accompanied by random 
jumps among the stable states induced by thermal or 
quantum noise. These large amplitude fluctuations have 
small probability away from the bifurcation points and 
the parametric threshold, but become significant in the 
vicinity of these critical points (cf. Ref. dl and references 
therein). These effects are out of the scope of this paper, 
here we restrict to small quantum fluctuations around 
well defined classical states outside of the critical regions, 
both below and above the parametric threshold. 

The bifurcation of the DufRng oscillator response is 
employed in Josephson bifurc ation amplifiers (JBA) for 
dispersive qubit readoutPSl^U, This method also applies 
to the parametric regime below the threshold (Josephson 
parametric bifurcation amplifier, JPBA). The novel fea- 
ture here is the possibility to measure amplitude of the 
amplified probing tone, which exhibits strong dispersion 
with respect to the detuning near the threshold, and can 
be advantageous for high fidelity qubit readout. 

The parametric radiation above the threshold offers 
yet another strategy for the qubit readout based on the 
significant contrast between the strengths of the output 
radiation above the threshold and the amplified noise be- 
low the threshold. 
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The paper is organized as follows. Sections p) III 



and IV are devoted to the development of the theoret- 
ical framework for describing parametric resonance in a 
high quality tunable cavity. In Sec. |V] we consider the 
nonlinear cavity response to a classical input signal be- 
low and above the parametric threshold, and in Sec. |VI| 
apply the results for the analysis of parametric amplifi- 
cation and methods of dispersive qubit readout. Section 



VII is devoted to the analysis of quantum fiuctuations. 



II. CIRCUIT LAGRANGIAN 

The device we study is sketched in Fig. [T] Its main 
part is a tunable superconducting strip line cavity termi- 
nated with a SQUUP^^. The cavity is weakly coupled 
to a transmission line that feeds an external microwave 
signal in and provides means for probing the field inside 
the cavity. The cavity is a spatially extended system of 
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FIG. 1. Sketch of a tunable cavity device: the A/4-cavity 
is terminated by a dc SQUID at the right end, and is ca- 
pacitively coupled to a transmission line at the left end; the 
SQUID is flux biased (phase 2/) via inductive coupling to 
a flux line imposing a driving phase 2/cxt; <t>d = (t>{d,t) and 
(j)o — <?!>(+0, t) are the phase values at the right and left ends 
of the cavity, respectively. An incidental signal fed in from 
the transmission line is reflected, separated from the input, 
and then analyzed. 



length d with inductance Lq and capacitance Co per unit 
length, and the cavity state is characterized by the su- 
percondu cting p hase field (j){x, t). We use the Lagrangian 
formalisrrP^IIll] describe the nonstationary dynamics 
of 0(x, t). 

The Lagrangian of the entire device consists of the sum 
of the Lagrangians of the cavity, transmission line, and 
the coupling, 



C[(t>] = £cav + CtL + 



(1) 



The Lagrangian of the cavity in its turn consists of the 
Lagrangian of the bare cavity £cav, and the Lagrangian 
of the SQUID CsiM, 
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dx (</)2 - w^^'^) (2) 
l + 2Ej cos f{t) cos (f>d 



Here v = I/v^IqCo is the field propagation velocity, 
(t^dit) = 4'{d,t) is the boundary value of the cavity field 
at the SQUID, and f{t) is the phase across the SQUID 
controlled by external magnetic flux, see Fig. [l] The 
SQUID is assumed symmetric for simplicity, with two 
identical Josephson junctions, each having a Josephson 
energy Ej and a capacitance Cj. The phase f(t) ap- 
pears in Eq. ([2| as an external time-dependent parame- 
ter that is able to excite parametric resonance. In fact 
it is a dynamical variable that describes, together with 
the variable (pd, the dynamics of two coupled Josephson 
oscillators of the SQUID driven by the external electro- 
magnetic field /cxt (t) ■ In Appendix [A] we show that in 
the limit of small 0^ ^ 1 the /-oscillator decouples from 
the (/)d-oscillator. Moreover, for experimentally relevant 
circuit parameters, the /-oscillator follows the drive field 
adiabatically because the resonance frequency of the /- 
oscillator is large compared to a typical resonance fre- 
quency of the cavity. A detailed derivation of Eq. ^ 
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and the connection of f{t) to the external field /oxt(i) is 
provided in Appendix [A] 

We assume here that the controlling field f{t) is com- 
posed of a constant biasing part F and a small harmonic 
oscillation with amplitude (5/^1, 



f{t) ^ F + Sfcosnt. 



(3) 



It is worth mentioning that the constraint 0^ ^ 1 is 
essential, otherwise the two Josephson oscillators become 
coupled and exhibit complex, even chaotic behavior un- 
der external driv^^H^ 

Proceeding to the other components of the device, we 
suppose the transmission line to have the same charac- 
teristic parameters Cq and Lq as the cavity. 
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The capacitive coupling is described with the Lagrangian 
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where (f)Q — <p{+Q,t) and cf)TL,o — 0tl(— 0,t) are the field 
values at the different sides of the coupling capacitor Cc- 



III. PARAMETRIC DYNAMICS OF CLOSED 
CAVITY 

We first consider the cavity decoupled from the input 
line, Cc = 0. The goal will be to identify the cavity 
frequency spectrum and investigate the parametric reso- 
nance. 



A. Cavity modes 

The Lagrangian £cav, Eq- explicitly contains two 
dynamical variables, the phase field (j){x,t), and its 
boundary value (j)d{t)- Variation of the associated action 
with respect to (j){x, t) leads to the wave equation. 



v'^" = , 



(6) 



supplemented by the boundary condition (f)^ — Q aX the 
open end of the cavity. Variation with respect to the 
boundary value 4>d{t) yields the boundary condition. 



Er 



2Ej cos f{t) sin (j)^ + 



0, (7) 



where Ec = (2e)V(2Cj) and El,^^^ = {h/2ef{l/LQd). 

Under static biasing, 5f — 0, the linearized boundary 
condition of Eq. ([t]) determines the set of cavity eigen 
mode^. 
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FIG. 2. Cavity spectrum: (a) graphical solution of Eq. Q 
for 1/7 = 25, and Cj = (solid) and C = 2Cj/{Cod) = 
0.05 (dashed); (b) cavity spectrum knd vs. I/7 according to 
Eq. ([9|, the vertical line indicates the value I/7 = 25 used in 
(a). 



The frequency spectrum a;„ is non-equidistant, and can 
be tuned by varying the bias F. 

Although the first term at the rhs of Eq. ([9| can 
in principle be tuned to zero, at = 7''/2, in prac- 
tice it dominates over the second term, at least for the 
lowest cavity modes, by virtue of the large parameter 
LUj/ujn ^ 1, where ljj — \j2EjEc jh is the Joseph- 
son plasma frequency. Indeed, given typical experimen- 
tal values, Ejjh w 4500 GHz and Ecjh - 10 GHz, the 
plasma frequency is wj k, 300 GHz, while the cavity fun- 
damental frequency is ~ 40 GHz, i.e., by one order 
of magnitude smaller (for typical cavity parameter^SHHl 
io ~ 4 • 10-^ H/m, Co - 2 • 10-^° F/m, and d « A/4). 

Furthermore, the cavity inductive energy is typically 
small, El^cscv/^ ~ 400 GHz, compared to the Josephson 
energy 2Ej. Taking advantage of this relation, and ne- 
glecting the capacitive term in Eq. (|9|, we get the ap- 
proximate solutions 
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(10) 
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The solutions of the spectral equation ^ are graphically 
illustrated in Fig. [2]ja), while Fig. [2]^b) shows the cavity 
spectrum as a function of the parameter I/7. 



B. Cavity Hamiltonian 

The Lagrangian formalism is sufficient for analyzing 
the classical parametric resonance. To describe the quan- 
tum dynamics the Hamiltonian approach is more con- 
venient. We derive the cavity Hamiltonian by expand- 
ing the cavity field over the complete set of cavity eigen 
modes, 



2c / 2 

0(a;,t) = ^^(?«(t)cosfc„a;, (12) 



where qnif) are time-dependent coefficients, and fc„ obey 
Eq. Using expansion (12) and noticing that the set 
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of functions cos fc„a; is non-orthogonal, we present the 
Lagrangian ([2| after some algebra in the form, 



Here the "masses" of the mode oscillators are given by 
the expressions. 



sin 2knd AC, 
2kad ^ dCo 



+ 37^cos2fc„d, Mo«l, (14) 



and 



V{qn, t) = -2Ej [cos /(t) cos c^d + COS (15) 

is a nonstationary nonlinear potential that mixes the 
eigen modes (see Appendix [b] for details of the deriva- 
tion) . 

It is convenient to absorb the factors Af„ and w„ into 
the rescaled coordinate. 



n^n Qn 



(16) 



and redefine the mode expansion in Eq. (12 1 accord- 



ingly. Then introducing the conjugated momenta, p„ = 
d£/dqn = Qn/i^n, wc arrive at the cavity Hamiltonian, 



C. Resonance approximation 

For small pumping amplitudes and weak non-linearity, 
the potential V{qn,t) in Eq. (17) could be considered 
perturbatively. However, the perturbative approach does 
not apply to the case of parametric resonance, when the 
pumping frequency matches an algebraic sum of the cav- 
ity eigen frequencies, 17 a;„ ± uj„i. In this case the 
corresponding cavity modes are strongly mixed and un- 
dergo complex time evolution. A particular case is the 
degenerate parametric resonance for m = n. In this pa- 
per we consider for certainty the degenerate paramet- 
ric resonance of the fundamental mode, « 2tJo- The 
method outlined below is straightforwardly extended to 
a non-degenerate parametric resonance. 

First we perform a canonical transformation corre- 
sponding to a transition to the rotating frame with fre- 
quency fl/2. This is conveniently done in terms of a 
complex variable. 



an = {qn + ipn)/v2h 



(18) 



for which the transformation reads a„(t) = e *^*/^A„(i). 
The equations of motion for the amplitudes An{t) read. 

An = -KcOn - n/2)An - ^Jl^^^^ . (IQ) 

^/2/l dqn 



At this point we take advantage of small values of the 
pumping amplitude, 5/^1, and the field amplitude, 
^ 1, and expand the potential Vijimt) in powers 
of these small parameters, keeping only the first non- 
vanishing terms, 



V{qn, t)Ki- {EjSf sin F cos nt) 



E 
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' cosF4. (20) 



Close to the resonance, $7/2 — cjq = (5 <C w„, a;„ — a;,„, 
the variable Ao depends slowly on time while all the other 
variables contain rapid time oscillations. After averaging 
over these oscillations we arrive at the shortened equation 
of motion for Aq (we skip the mode index below). 



A -id A- ieA* - ia\A\^A = , 
with the parameters 

SfuQ tan F cos^ kod 

"2 knd 
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(21) 



(22) 



(23) 



When applying the canonical transformation to coordi- 
nate and momentum, A = {Q + iP) / \/2h, and averaging 
over fast oscillations, the cavity Hamiltonian is cast into 
the form, 



a 
8h 



.(24) 



This Hamiltonian corresponds to the metapotential of 
the parametric lumped element oscillator^, i.e. the de- 
generate parametric resonance in the cavity is mapped 
on the one in a lumped element oscillator. The map- 



ping is defined by Eqs. (22) and (23), where the effective 



pump strength e, and the nonlinearity coefficient a are 
expressed through generic cavity parameters. 

According to the experimental values discussed in 
Sec. HI A parameter 7 in Eq. (10) is estimated as 7 ~ 
4 • lO"'^. For such a small value of 7, the parameters e 
and a are approximated, using the spectral equation (|9| , 

e/wo w 7(5/ tan F/2 < 5f , (25) 
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< 10- 



(26) 



where Zq = y/Lo/Co is the cavity impedance and Rq = 
h/2e^ is the quantum resistance. 

It follows from these estimates that the effective pump 
strength e is substantially reduced compared to the am- 
plitude of the phase modulation in the SQUID, and the 
effective nonlinearity of the cavity oscillator is signifi- 
cantly smaller than the underlying bare nonlinearity of 
the SQUID oscillator (a = wj/6 for the Josephson po- 
tential) . These remarkable properties result from the fact 
that the cavity is almost shortcut to the ground at the 
edge X — dhy virtue of large Josephson energy in Eq. ([?]) 
(7 <C 1), hence the boundary value of the field amplitude 
(bri is small. 
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The small values of the effective oscillator parameters 
are essential for the validity of the resonance approxi- 
mation. The latter requires the evolution of A to take 
place on a time scale much larger than the period of the 
cavity fundamental mode, 1/ujq, over which the initial 
Hamiltonian is averaged. 

It is instructive to express the constraints earlier im- 
posed on the phases, (f>d,Sf ^ 1, in terms of the ampli- 
tude A and the pump strength e. 



\A\ < sJrJz^, e < 7Wo . 



(27) 



or equivalently, oLAp, e ^ ^ujq <^ ujq. In other words, 
the constraints ( |27[ ) are more stringent than the ones re- 
quired for the resonance approximation, a|v4p, e ^ ujq. 
On the other hand, these constraints provide sufficient 
room for the pumping strength to be increased above 
the parametric threshold beyond the resonance width T 
(see Eq. (351 in the next section), a|Ap, e ~ T <C 7^9, 



for a high quality cavity. 

In most of our calculations we restrict to the lowest 
order ^/-dependence in Eq. ([20]), however, in some cases 
it is useful to keep higher order terms. In particular, the 
second order term cx Sf^ will introduce, after averaging 
over time, a nonlinear shift of the resonator frequency, 
proportional to e^. This shift is evaluated in Eq. (All I in 
Appendix [Aj and in terms of the effective pump strength 
it reads, 



ujo (e) - Wo 



Wo 



juiq tan^ F 



(28) 



This shift could be used in practice for evaluating the 
actual magnitude of the pump power acting upon the 
SQUID, which is usually not known. Also, it causes 
quenching of the parametric instability at large pump 
strength, as will be shown in Sec. |V A| 



IV. CAVITY COUPLED TO TRANSMISSION 
LINE 

The parametric effect in the closed cavity is an ide- 
alization. The connection to the external transmission 
line gives rise to the qualitatively important new fea- 
tures: firstly, the cavity field is allowed to leak out of the 
cavity, giving rise to the cavity damping, and secondly, 
an external electromagnetic signal can be fed into the 
cavity and amplified. Our aim in this section will be to 
include these features into Eq. (21), and derive the rela- 



tion between the input and output fields, thus preparing 
the framework for the further investigation of parametric 
amplification. Ou r derivation closely follows the input- 
output theorjJsElI^ ^ggg also illuminative derivations in 
Refs.[5]and|^. 

Aiming at the analysis of the quantum dynamics of 
the open cavity, we describe the field in the transmission 



line in terms of spatial modes, similar to Eq. ( 12 1 for the 



cavity, 

4>TL{x,t) 
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dk 



qk{t) cos kx, (29) 



with Wfc = vk. Opening of the cavity invokes also an 
additional set of modes oc sin kx, however, in the weak 
coupling limit these modes do not contribute in the main 
approximation and are neglected here. 

Focusing on the effect of cavity damping at weak cou- 
pling, we will only keep the cross term in the coupling 
Lagrangian ([s]). 



C, 



c<PO(PTL,0 ) 



(30) 



and neglect the quadratic terms, thus neglecting small 
corrections to the kinetic energies. With this simplifica- 
tion, and retaining only the fundamental mode field in 
the cavity Lagrangian, we write the total Lagrangian in 
the form. 
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1 
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(31) 



The corresponding Hamiltonian reads, to first order of 
the weak couphng (Cc <C Cod), 



-V{q,t) 



dkujk {pI + ql) 
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(32) 



Repeating the derivation of the previous section we 
derive coupled equations of motion for the cavity ampli- 
tude a and the spectral amplitudes of the transmission 
line, flfe = {qk + ipk)/V2h, 



id = LOoa + 
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(34) 



Here we introduced the cavity damping rate. 



To = ujo 
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(35) 



Near the parametric resonance the equations of motion 
for the slow variables, A{t) — e'^*/^a(t), and Ak(t) = 
e'^*/^afc(t), take the form, after averaging over rapid time 
oscillations. 



iA+{S + a\A\^)A + eA* = 
lAk-SkAk = ^I^^A. 



irko 



irko 



dky/u^Ak (36) 



(37) 
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with Sk = ujk — VL/2. 

We ehminate the transmission hne modes from 



Eq. (|36|), invoking the solutions of Eq. (37) 



Ak{t) = Afc(to)e-'*''=(*-*") (38) 



Trfco 



with initial conditions Ak{t{j) at time to < and substi- 
tuting it into Eq. ( 36 1 . Within the resonance approxima- 
tion, the factor ^^w^ in the integrand is to be replaced 
with y^wq, and the integration over the wave vector k to 



be extended to the entire axis. After making these ap- 
proximations we arrive at the Langevin equation for the 
cavity amplitude. 



iA + 5A + eA* + a\A\^A + iEo^ 



(39) 



with the input flux amplitude 
1 



B{t) 



d4A,(io)e-'*^(*-*°). (40) 



The amplitude B{t) is associated, as shown in Appendix 
[C| with the incident (right-going) wave in the transmis- 
sion line, B{t — x/v), taken at the boundary a; = 0. 

The solution of Eq. ( 37 ) can be equivalently expressed 
in terms of the amplitude at a future time, Ak{ti), ti > 
t, which defines the output flux amplitude C{t) via a 
relation similar to Eq. (40) with ti substituting for tg. 



This output amplitude is associated with the reflected 
(left-going) wave in the transmission line, C{t + x/v), 
taken at a; = (Appendix C]). The relation between the 
output and input amplitudes reads, 



C'it) = B{t) - iy/2roA{t) 



(41) 



The parametric pumping couples the cavity field am- 
plitude A and its complex conjugate, and it is convenient 
to rewrite Eq. (39) in the matrix form. 



(42) 



where 



A 

c 



C + ir e 
e C-ir 

^a\A\\ 



(43) 



The conservative part of the dynamics in Eqs. (39 1, 



( 42 ) is determined by the effective Hamiltonian 
H{Q,P) ^ H,,, + 2^o\B\{cos eBQ + sin Ob P), (44) 



with HcuviQ, P) from Eq. (24), and 6b being the phase 
shift between the input amplitude B — |_B|e'^^ and the 
pump. 

Besides the damping Tq associated with the opening of 
the cavity, there might also be internal losses in the cav- 
ity, e.g. caused by the cavity resistance. A way to account 



for these losses is a model with a fictitious transmission 
line coupled to the cavity, that acts as a scattering chan- 
nel with a noisy input amplitude Bj^it) and an associated 
damping rate Ffl . This would lead to an enhanced damp- 
ing rate, F = Fq -I- F^, at the Ihs of Eq. (39 1, and also 



introduce an additional input term, \/2TjiBii(t), at the 
rhs of this equation. 

The damping effect results in the broadening of the res- 
onance, and if the resonance becomes sufficiently broad, 
higher cavity modes might also be excited, despite the 
non-equidistant property of the cavity spectrum. In this 
case, the isolated mode dynamics of Eq. ( 39 ) would be 



replaced by a more complex dynamics of parametrically 
excited coupled modes. To ensure the validity of the 
single-mode approximation, the condition \iOQ ± F/2 -|- 
ri — Wi| > Fi/2 must be met, where Fi is the reso- 
nance width of the first cavity mode. For ft — 2ujq, 
and the cavity spectrum given by Eq. ^ and param- 
eters of Sec. |III A[ the anharmonicity is of the order, 
ujQ + il/2 — u!i = 'Sluq — oji ~ This implies that 

the cavity quality factor Q — ujq/T should not be less 
than Q > 10"^ . This corresponds to a small coupling ca- 
pacitance in Eq. (35 ), Cc/CqCL ~ 10~^, assuming that the 



internal losses are not dominant, Tr < Tq. 



V. CLASSICAL CAVITY RESPONSE AND 
RADIATION 

In this section we analyze the cavity response to a 
noiseless classical input signal. We consider harmonic 



,-iAt 



where 



inputs, which have the form B{t) = Bi 
A = Lo — rt/2 is the detuning of the input signal from 
the half frequency of the pump. For the input frequency 
A = the cavity response is stationary, and it can be 
fully analyzed in the nonlinear regime. For detuned in- 
puts, we restrict to small input amplitudes; at large am- 
plitudes the nonlinear response becomes complex and ex- 
hibits a transition to a chaotic regime. 



A. Parametric resonance in absence of input signal 

We start with the analysis of the intrinsic parametric 
resonance in the cavity in the absence of input signals, 
B{t) = 0. Due to the damping, any initial cavity state 
evolves towards one of the steady states that define the 
picture of the parametric resonance. These steady states 
depend crucially on the pump strength e, and also on 
the detuning of the pump frequency from the cavity res- 
onance, (5 = f2/2 — cjQ. 

If e < F, only the trivial steady state, ^ = 0, exists for 
all values of the detuning i5. If e > F, the trivial state 
turns unstable within the interval \5\ < — F^, and 
instead two non-trivial stable steady states, A — |^|e'^-*, 
emerge at the threshold 5 = Ve^ — T^j and persist for all 
5 < Ve^ — F^, see Fig. K5la). These states have identical 
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FIG. 3. (a) Amplitudes of cavity steady states vs. de- 
tuning S for e = \/2r, stable (solid) and unstable (dotted); 
dashed vertical lines separate mono-, bi- and tristable re- 
gions; insets show phase portraits of corresponding regions, 
(b) Boundary of parametric instability without (solid) and 
with (dotted) account of nonlinear frequency shift, Eq. (28 1, 
for a = r/100 and ^ = 1/10, cf. Eq. (|47|; yellow region cor- 
responds to bistable high-amplitude state, blue region indi- 
cates coexistence of stable high-amplitude and zero-amplitude 
states, (a = T/lOO). 



amplitudes, 



(45) 



and are 7r-shifted in phase, with sin(20^) =T/e. 

In the further red detuned region, S < —y/e^ — V^, 
the trivial steady state solution, A — 0, becomes stable 
again, such that the three stable states coexist there. 
Simultaneously, two new unstable states emerge having 
the same amplitude, \A\'^ — (—(5 — Ve^ _ ^q,^ 

In the limit of e = 0, F = (the undamped Duffing os- 
cillator), the nontrivial stable and unstable states merge, 
forming a manifold of marginally stable states with in- 
definite phase 9a and amplitude |Ap = —S/a. 

The steady states of the damped cavity at e > F origi- 
nate from the fixed points of the cavity Hamiltonian -ffcavj 
Eq. (24), which are illustrated in the insets of Fig. [sj^a) 
for the mono-, bi-, and tristable regions. The damping 
F introduces the threshold for the emerging nontrivial 
states, and shifts the positions of the steady states in 
phase space away from the fixed points. 

The pump parameters where new steady states occur 
are determined by the stability properties of the underly- 
ing linear system, characterized by the matrix A{a — 0). 
Its determinant, D = causes divergence at the 

parametric instability threshold, |^| — \J — T^, where 
the fixed point A = turns unstable. In a linear system 
this would lead to exponentially growing solutions in the 
parameter regime e > F and \b\ < \J ~ T^, with a rate 




-100 



FIG. 4. Relative areas of attractor basins of coexisting high- 
amplitude (lower curve) and zero-amplitude (upper curve) 
stable states vs. detuning Insets show the basins of attrac- 
tion of high-amplitude (red and green) and zero-amplitude 
(blue) states, for (5/r = -30 (left) and (J/F = -4 (right) 
(e = V2r, a = F/lOO). 



A = — F + \J — . In the nonlinear system this global 
instability is lifted by the bifurcation of the fixed point 
A = into the two new stable steady states. 

The cavity field, as it leaks into the transmission line, 
generates an outgoing field with the amplitude C accord- 
ing to Eq. ([4T|). For the steady state, Eq. ([45]), the flux 



radiated into the transmission line amounts to 
2Fo 



\CY 



-S + - F2 



(46) 



The nonlinear effect of the cavity resonance shift in- 
duced by the pump, mentioned in Sec. IIIC Eq. (28 1, 
leads to the quenching of the parametric instability at 
strong pumping as observed in experimenlPS-. By tak- 
ing into account this shift, the actual pump detuning be- 
comes n/2 — ujo{e) = S— {ujQ (e) — wq), and the parametric 
instability condition modifies accordingly. 



S<^€ 



Y2 



F '' 



(47) 



where /? — T/uJojtan^ F <C 1. The modified boundary 
of parametric instability in the {S, e)-plane is depicted in 
Fig. [Sj^b): the instability region is bounded by the maxi- 
mum blue detuning, 6ma.x = r/4/3, and it is also bounded 
by a maximum pump strength at given detuning, e.g. 
Emax = T/p at 5 = 0. 

In the experiments^, all the described states of the 
parametrically pumped cavity have been observed: the 
subthreshold monostable regime at blue detuning, as well 
as the above-threshold bistable and tristable regimes at 
red detuning. The visibility of particular stable states in 
the multistable regime is defined by the probabilities of 
their occupation, which are determined by the relative 
areas of the respective basins of attraction, i.e. the phase 
space regions from which trajectories asymptotically ap- 
proach the respective state. Examples of the attractor 
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basins in the red-detuned region, S < — -v/e^ — F^, are 
shown in the insets of Fig. |4] where the blue basin be- 
longs to the zero-amplitude state, and the red and green 
attractor basins are those of the high-amplitude states. 
The relative areas of the latter rapidly decrease and be- 
come very small in the far red-detuned region, as shown 
on the main panel in Fig. [4j implying that these states 
are much less populated. A similar conclusion is drawn 
from the calculation of the probability to escape from 
the high-amplitude stateJ^l, which is much larger than 
the one for the trivial state, A = 0, at far red detuning. 

These arguments explain why in the experimenli^SJ the 
boundary of parametric resonance is washed out at red 
detuning, in contrast to the sharp boundary at blue de- 
tuning, which is determined by the threshold for the non- 
trivial steady states. 



B. Driven Duffing cavity (e — 0) 

Now we turn to the discussion of the cavity response 
to a weak signal with zero detuning, A = 0, and complex 
amplitude B = \B\e'^^. 

It is instructive to first review the response of the 
driven Dufhng oscillatoJ^', which corresponds to the limit 
e = in Eq. (42). In this case the detuning S refers 



to the deviation of the input frequency from the cavity 
resonance. The cavity response is given by the equation. 



B. 



(48) 



The maximum response is achieved at C = 0, along the 
tilted line |Ap((5) = — (5/a, and amounts to |AmaxP = 
|^|2(( = 0) = 2Fo|B|Vr2, independent of a. As a 
consequence of the tilted resonance line, the cavity re- 
sponse can display bistability, with two coexisting stable 
states, as shown in Fig. |5ja). The bistability emerges 
above the critical value of the driving amplitude, \Bc\^ = 
4F3/(3V3aro), and at the detunings, 5 < 5c = -V^r. 
The bistability region is confined by the bifurcation lines. 



\B+\ = 



27aTQ 



-1 - 



9r2 

1^ 



± 1- 



3r2 

1^ 



3/2' 



(49) 



forming a wedge in the ((S-|i3p) plane, as illustrated in 
Fig. [6] with black lines. 

An ideal Duffing cavity fully reflects the input signal, 
so the amplitude of the output, C = |C|e'^"^, carries no 
information about the resonance, |C| = \B\. Such infor- 
mation is only available for a lossy cavity, where 



\cl 

IBP 



1 



4rnr 



0^ R 



(50) 



On the other hand, the phase 9c' of the output sig- 
nal is sensitive to the position of the resonance. This 
is the working principle of the Josephson bifurcation 
amplifier^^I!^ where the variation of 9c under sweeping 




FIG. 5. Steady state cavity response \A\'^ according to 
Eq. I\5l\ vs. pump detuning S for different values of the pump 
strength e, below threshold, e/T — 0.9,0.93 (a,b), and above 
threshold, e — 1.1, 1.5 (c,d). Solid and dashed lines mark sta- 
ble and instable states, respectively. (j-Bp — 2r, 6b ~ 7r/2, 
a = r/100, Tr^O). 
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FIG. 6. Bistability regions of the cavity response, bounded 
by bifurcation lines |i3±(5)p, for different subthreshold values 
of the pump strength, e/F = 0,0.4,0.7, 1.0 {Ob = tv/2, a = 
r/100, Fn = 0). 



the input power through the bistability region is ex- 
ploited for the qubit readout. 



C. Driven parametric cavity 

Switching on the parametric pumping, e > 0, qualita- 
tively changes the cavity response. Now the amplitude 
of the cavity field is determined by the equation. 



1^1 



F^ + £2 _ 2e[C cos 261b F sin 26*5]) , 
D = det{A) = -f- F^ - £2 . (51) 
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In the subthreshold regime e < F, the cavity response 
remains quaUtatively similar to the Duffing oscillator, see 
Figs.[5ja)-(b). The role of the parametric pumping in this 
regime is to effectively reduce the damping term, — >■ 
— e^. This makes the resonance more narrow and, at 
the same time, strongly increases the cavity amplitude 
along the tilted resonance line C = 0. Another important 
feature is an explicit dependence of the cavity field on the 
phase shift 9b oi the input with respect to the parametric 
pump. 

The maximum value of the cavity field is. 
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,C = o) = 2ro 



(r2 



(52) 



Similar to the Duffing limit, this value is independent of 
the nonlinearity coefficient a. The maximum response di- 
verges at e = r, which can be compared to the resonance 
catastrophe of a linear parametric oscillator. While in 
the linear case the divergence occurs at (5 = iVe^ +r2, 
the nonlinearity here shifts the divergence towards an in- 
finite red detuning. 

As a consequence of the resonance narrowing, the crit- 
ical bifurcation point moves towards the origin, |i?cP — 
dc = when e — )• F, as illustrated in Fig. [6| 

Above the threshold, e > F, the resonance splits into 
two branches, as shown in Figs. [5f^c)-(d), each branch 
consisting of two non-degenerate steady states, one pair 
being stable and the other unstable. These states origi- 
nate from the degenerate nontrivial states in the absence 
of an input signal, cf. Fig. [3j the degeneracy being now 
lifted by the input. The distance between the branches 
increases with e. 

The scattering by the parametrically pumped cavity is 
always inelastic, in contrast to the Duffing cavity, and the 
output signal in general differs significantly from the in- 
put signal, not only in phase but also in the absolute 
value, \C\ 7^ \B\. Using the input-output relation in 
Eq. ([ilj), and the steady state solution \A\'^ in Eq. ([42|), 



the output amplitude can be expressed as a function of 
the input amplitude. 



C 

C* 



V 



B 

B* 



que 



"I 



-iq 



with the parameters 
2eFo 



q = 



D ' 
arctan 



qR = 



Vi + q 

4FoFfl 

D 
'2FoC 



1? 



qne 



-"I 



(53) 

(54) 
(55) 



The relation in Eq. ( 53 1 maps the points of the unit circle, 
B = e'^^, onto the phase-dependent curve, C{9b), and 
determines the phase-dependent gain 

\ci9Br 



= I + 2g2 



\B\ 



qR 



- qRs'm{29B +?/) .(56) 
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FIG. 7. Anisotropy of the cavity output field in the complex 
C-plane, in the subthreshold regime (e/F = 0.2,0.7 [blue, 
red]) and above the threshold (e/F — 1.4 [yellow]); insets 
show the dispersion of the gain with the input phase 6b below 



threshold (left) and above (right). {S 
F/lOO, Fr = 0). 



0, 



2r, a 



The 6'B-dependence of the gain ICp/jSp and the quadra- 
tures of C = (A -f iY) /2 are illustrated in Fig. [t] In the 
monostable (subthreshold) regime the output amplitude 
C is amplified (G > 1) or deamplified (G < 1) depending 
on the input phase. For e < F the points C{9b) form 
a strongly elongated curve in phase space, centered at 
(0, 0). In the quasilinear regime, where the parameters q 
and qn in Eq. (56) are approximately independent of 0^, 



this curve approaches an ellipse with the half axes 



(57) 



giving the maximum / minimum gain factor along those 
quadratures. For negligible internal losses, qn/q^ ^ I, 
the amplified and deamplified quadratures are related ac- 
cording to Ginin = 1/G,nax- 

In the limit e — the gain factors become equal and 
reduce to the reflection coefficient of the Duffing oscilla- 
tor, Eq. (l50). 



In the bistable regime above the threshold the corre- 
sponding output amplitudes C{9b) are mapped on two 
distinct closed curves in phase space, with a 7r-phasc shift 
between them, as shown in Fig. [7] The offset from the 
origin is due to the parametric radiation generated by 
the cavity. 



D. Response to detuned signal 

The cavity response has a simple stationary form only 
when the frequency lo of the input signal strictly matches 
the half- frequency of the pump, 11/2. If the input is time- 
dependent in the rotating frame, e.g. B{t) ~ Ber^^*' 
with A = w — r2/2, the combination of the time-periodic 
force with the nonlinearity leads to the formation of a re- 
gion in phase space where the cavity amplitude A evolves 
chaotically, as illustrated in Fig.[8]for the bistable regime 
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J± u 



FIG. 8. Phase-space representation of cavity amplitude A = 
(Q + iP)/y/2h driven by detuned input amplitude B{t) = 
Be-'^\ The trajectories {Q{t),P{t)) are evaluated in the 
conservative limit, neglecting the damping term on the Ihs of 
Eq. ([39|, for (a) A = 0, (b) A = T/IO, and (c) A = T. In 
(b)-(c) trajectories are represented stroboscopically, at times 
t = 27m/A, {n = 0,1,2,...). (e = V^T, 5 = 0, l^]^ = T, 
9b = n/2, a = T/lOO, Tr = 0). 
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FIG. 9. Linear gain of detuned input signal; (a) G(A) 
for small detuning S — l.SSF (single resonance), and large 
detuning 5 — 7T (split resonance), at fixed pump strength 
e/r — 2; (b) G(A) for fixed detuning 5 = SF and increas- 
ing pump strengths, e/F = 3 . . . -^/l -|- S'^/F^ (from bottom to 
top). {Tr = 0). 



above the parametric threshold. With increasing input 
amplitude and detuning a chaotic layer forms around the 
instable fixed points of the Hamiltonian ( 44 1 , and its area 
grows with \B\ and |A|, see Fig. [sj However, as long 
as the stable fixed points persist in the presence of the 
time-dependent drive, the amplitude of the damped cav- 
ity evolves into a time-periodic limit cycle around them, 
and then the time-average of A{t) gives only small cor- 
rections to the stationary result. 

In this section we evaluate the response of an ideal 
cavity (F/? = 0) to a detuned signal in the monostable 
regime, e < V S'^ +T'^. We restrict to a linear response 
assuming a|^P <^ + S"^ — e"^. 

Suppose the input signal in Eq. ( 39 ) consists of two 
conjugated harmonics, B{t) = B(A)e""^* + B(-A)e'^* 
(signal and idler in the terminology of non-degenerate 
parametric amplification). Then the output field, as well 
as the field in the cavity, will also consist of the combi- 
nation of the same harmonics. The output amplitudes 
are related to the input via the equation generalizing 
Eqs. (pi)- 154 



/C(A) 
C*(-A) 



where 



V(A) 



1 



«ii(A) 

^(A) = (r-iA)2 



«ii(A) 
D{A) V «i*2 ^^n(A) 
(S - ir)2 - A2 '2 



Vl2 



e , ui2 = 2ire , 



(58) 



(59) 



(60) 



The coupling between the conjugated harmonics is a fin- 
gerprint of parametric amplification: an input at fre- 
quency A generates outputs at frequencies A and —A, 
and conversely an output at frequency A consists of the 
contributions of inputs at frequencies A and — A. In par- 
ticular, for B{-A) = 0, Eq. ^ yields. 



(61) 
(62) 



C{A) = v/l + k(A)|2ei''(^)B(A) , 

C{-A)=iq*{A)B*iA), 

q(A) = 2ero/D(A). 



Amplification of the detuned signal is character- 
ized by two gain factors, direct gain Gi(A) — 
|C(A)|2/|B(A)|2 = l + |g(A)|2, and interconversion gain, 
G2(A) = |C(-A)|V|B(A)|2 |g(A)|2. These two gains 
are fundamentally related, G2 = Gi — 1, which is the 
consequence of the fundamental property of the matrix 
elements in Eq. ([59|, \vii\'^ - |wi2p = \D{A)\^. For the 
quantum fields, this property guarantees the unitary re- 
lation between the input and output quantum states (see 
later in Sec. 



VII) 



The amplification of detuned signals possesses another 
interesting property - the appearance of resonance fea- 
tures, as illustrated in Fig.|9][a). The resonance structure 
of the gain is determined by the determinant |_D(A)p, 
Eq. (62). It has a single minimum, at A = 0, within 



the interval of relatively small detuning, F^ 

-f r^, and the gain factor G(A) is accordingly single 
peaked at A = 0. However, at larger detunings. 



(63) 



two resonance peaks emerge, situated symmetrically with 
respect to A = at 



A = ±y/5^ - £2 - F2 



(64) 



The origin of these resonances can be understood from 
the behavior of the response function of a conventional 
damped linear oscillator, x('^) = i^o — — iuT)^^. At 
small damping, F <C wg, the resonance is close to the 
eigen frequency ujq, uj — ujq — /2. With increas- 
ing damping the resonance is pulled towards the zero 
frequency, and stays at the zero frequency as soon as 
F > V2ujq. Similarly, the resonances in the response 
of the linearized parametric oscillator, Eq. (39 1 with 
a = 0, are at small F close to the oscillator eigen fre- 

^2 



e^, as in Eq. (64), but are pulled to- 



quencies, ±VS^ 
wards A = with increasing F, and eventually merge 
when F > V(52 -e^, Eq. ([63|. 
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VI. AMPLIFICATION AND QUBIT READOUT 

In this section we discuss the apphcation of the para- 
metrically pumped cavity for signal amplification, and 
for dispersive qubit readout. 

In what follows we shall neglect internal losses in the 
cavity and assume Fq = F. 



A. Amplification 

The amplification characteristics of the nonlinear para- 
metric cavity depend on many parameters: pump and in- 
put strengths and dctunings from the cavity resonance, 
relative phase shift, nonlinearity and damping, which 
makes the overall picture pretty intricate. 

The output power |Cp as a function of the input power 
for on-resonance input, A = 0, is depicted in Fig. 10 
for various values of pump strengths and pump detun- 
ings. The major phenomenon here is the appearance of 
multistable regimes. The bistable regime establishes al- 
ready below the threshold, e < F, in the red detuning re- 
gion, (5 < 0, as shown on Fig.[Io]^b). Above the threshold, 
the mono-, bi-, and tristable regimes exist at different de- 
tunings, as shown on Fig. [lOjc). Moreover, in the latter 
regime, the output power does not approach zero value 
at = due to the effect of parametric radiation. 

For the amplification purpose the monostable regime in 
Fig. 10 a) is the most suitable. The output power in this 
regime depends monotonically on the input power, but 
exhibits pronounced nonlinearity with increasing pump 
strength at input power levels \B\'^ ^ F. 

The maximum differential gain is achieved at small 
input power, and for phase shift ds ~ — 7r/4. The gain 
is controlled by the quantity q in Eq. (57), and at large 



AeVV 



> 1. 



(65) 



The gain increases while approaching the threshold 



(cf. Fig. ^h)), G„ 



(4eF)V(52 + r2 



the quasilinear approximation, and then it is limited 
by the nonlinearity. Let us evaluate this upper bound 
for the gain at e = F and ^ = 0. In this case, 
Gmax = (2F/a|A|2)4. Extracting the amplitude \A\'^ 
from Eq. (51]), with 9b = -7r/4, 



a|A|2 «F(8a|B|Vr' 



,1/5 



we get 



G„,ax = (4F/a)^/^ (|B|VF 



-4/5 



(66) 



(67) 



In a similar way we can evaluate the absolute minimum 
of deamplification. This is achieved at 0b « 7r/4, where 
G,„in « l/(2g)2, and 



«F(2a|B|Vr 



2/T-2\l/3 



leading to the equation for minimum gain. 



(69) 



We note that the nonlinear deamplification is more ef- 
ficient than the amplification: the product of the maxi- 
mum and minimum nonlinear gains significantly deviates 
from unity, in contrast to the linear case, 



G„,axG„,i„ « ^a\B\yAT^ < 1. 



(70) 



With these results we conclude that the maximum am- 
plification (deamplification) efficiency is controlled by the 
parameter F/a, and therefore a relatively small nonlin- 
earity coefficient is required to achieve a large parametric 
effect. 

As we will see later, the same conclusion is also valid 
for the nonclassical properties of the fluctuations. 

At this point it is appropriate to estimate the output 
signal-to-noise ratio for parametric amplification, refer- 
ring to the results of the noise analysis in Sec. VHP 
According to Eqs. (101) and (105) the amplified noise 



increases in the vicinity of the threshold, however, the 
noise amplification is less efficient than the signal ampli- 
fication, giving the ratio (for the quasilinear limit). 



\cl 

„vac 



\B\ 



(71) 



(68) 



This ratio is large as soon as \B\'^ > (F — e)/8. 

Amplification of a detuned signal, A 7^ 0, has qualita- 
tively similar properties in the vicinity of the parametric 
threshold, — (5^ ^ T^- Here the gain factor G(A) has 
a quasi-Lorentzian shape, peaked at A = 0, as shown in 
Fig. [9j^b) , the maximum gain increases while approaching 
the parametric threshold, while the bandwidth shrinks to 
zero. 

However, the bandwidth can be considerably in- 
creased, maintaining rather high gain, by working away 
from the parametric threshold in the region where the 
gain peak splits, 6^ ~ -f F^, Eq. (63 1. Here a wide 
frequency plateau emerges around A = 0, see Fig. [of^b), 
where the gain factor is nearly constant over a frequency 
interval given by the distance between the resonances, 
A = ±V(52 - £2 _ r2. 



B. Bifurcation readout below threshold 

The bifurcation regime of the cavity nonlinear response 
in the absence of parametric pumping is employed for 
dispersive qubit readout using JBA^, for a review see 
Ref. [3T] and references therein. With this method, the 
phase shift of a reflected (or transmitted) probing signal 
is measured while ramping the signal amplitude. The 
result is sensitive to the detuning of the signal tone from 
the cavity resonance, which is pulled by the qubit by ±5^, 
depending on the qubit state. 
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FIG. 11. Average output power (|Cp) vs. input power \B\^ 
for e/r = 0.7 and 5 /V = -1.7,-0.7,-0.5,0.5,5.0. Inset: 
maximum gain factor G vs. and detuning 5\ black lines in- 
dicate the boundaries of the bistable region (wedge) cf. Fig. [6] 
white vertical lines indicate the parameter traces used in the 
main figu re. The average output power (|C|^) is obtained 



FIG. 10. Output power |Cj vs. input power |S| , below input 
threshold, e < F (a,b), and above threshold, e = 1.2r (c). (a) 
5/V = 0.5 and e/F = 0, 0.3, 0.7, 0.9, 1.0 (from bottom to top); 
(b) 5/T = -0.72 and e/F = 0,0.63,0.7,0.8,1.0 (from bottom 
to top); the dotted line refers to the Duffing limit, e = 0. (c) 
e/F = 1.2 and 5/V = -1.4,-0.5,1.0 (from top to bottom) 
For each of the parameters instable branches are indicated by 
dashed lines. {63 = 7I-/2, a = F/lOO). 



One may take advantage of the high parametric gain 
for probing a qubit state by measuring the amplitude of 
the output signal instead of the phase shift. The ampli- 
fied signal exhibits significant dispersion over the cavity- 
pump detuning thus providing high contrast for the qubit 
readout. 

The basis of the method can be understood from 
Fig. [TT] here the average output power is plotted against 
the input power for different detunings below the thresh- 
old, e ~ 0.7 r. The bist ability wedge for this pump 
strength is illustrated in the inset, compare also Fig. |6] 
The lowest three curves in Fig. 11 correspond to values 
of the detuning within the monostable regions, either to 
the right or to the left of the critical bifurcation point, as 
indicated by white cuts in the inset (in the latter case, 
5 = — l.TF, the ramped input signal should not cross 
the bifurcation line). The other two curves correspond 
to crossing through the bistability wedge or very close 
to the critical bifurcation point (here the average output 
power in the presence of classical noise is plotted, which 
then exhibits a gradual transition from the low- to the 
high-amplitude branch of the bifurcation curve). 

The output contrast is extremely sensitive to the de- 
tuning: it is up to factor of 10 for detunings differing 
by a linewidth F already at rather small input power, 
~ 10 F. In practice, a cavity frequency pull ex- 



from Eq. ( 39 1 in the presence of white Gaussian noise in the 
1 = 7r/2, 



F/lOO). 




5e-l-02 



0.9 - 




le+03 

1 








1 

2e-01 




0.7 - 

-] 


.5 


-1.0 -0.5 0.0 





5 1.0 1 

S/T 


FIG. 


12. 


Output power |Cp/F 
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detuning 5 and 
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strength e for B = 0. The black line separates bistable and 
tristable regimes, cf. Fig. [3] the white line spans between two 
qubit-state dependent, effective detunings 5q ± 5q {5q — 0.7F, 
e > ^F2 + 52, So > Ve2-F2-5,, a = F/lOO). Inset: noise 



photon number n^'^/F, Eq. (101 1, vs. 5, e as in main figure 



erted by the qubit may be of the orde^^ 5q ~ 10 MHz, 
i. e. comparable to the linewidth, F 10~^ci;o < 10 MHz. 

The output contrast can be further enhanced by in- 
creasing the pump strength towards the threshold. It 
is also possible to ramp the pump strength rather than 
input power. The possibility to operate with several pa- 
rameters gives room for further optimization. 
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C. Radiation readout above threshold 

An alternative strategy for the dispersive qubit read- 
out is provided by parametric radiation above threshold. 
This method, illustrated in Fig. 



12 



is based on the fact 
that in the absence of an input signal, \B\'^ = 0, the 
output signal is zero in the monostable region below the 
threshold (at blue detunings), S + Sq > \J — T^, while it 
is finite above the threshold, |<5— 5g| < \J ^ — T^, where it 
equals, |C|2 = (2V I a) (-(5 - (5,) + Ve^ -r2), according 
to Eq. (Ilel). 



The maximum contrast is achieved by choosing the 
pumping strength, e > .y/(5^ + F^, and the optimum bi- 
asing detuning, b ~ Ve^ — F^ — (5^ , as illustrated in 
Fig. 12 Such a choice guarantees that the blue shifted 



point, (5 -I- lies in the monostable region close to the 
threshold, while the red shifted point, 8 — 8^, lies in the 
bistable region and not in the tristable region where the 
trivial cavity state, |^| = 0, dominates. Then the output 
radiation power does not depend on e. 



(72) 



This value is to be compared to the noise value in the 
monostable region below the threshold. The amplified 



vacuum noise is given by Eq. ( 101 1 in Sec. VII and illus- 
trated in the inset of Fig. [12] 



Fe^ 



r2 + ^2 „ £2 



B = 0. 



(73) 



Since the noise diverges at the threshold, the point 5 + 
is to be chosen not too close to the threshold. It is 
sufficient to depart from the threshold by ~ F to have 



the noise level. 



F. Then for 6„ ~ F, the radiation 



to noise contrast becomes. 



4F 



> 100 . 



(74) 



VII. 



QUANTUM FLUCTUATIONS OF CAVITY 
FIELD AND EMITTED FIELD 



So far we discussed the classical regime of paramet- 
ric resonance in the tunable cavity. In this section, we 
extend the formalism to the quantum regime, and in- 
vestigate the quantum properties of the field inside the 
cavity, and of the output field. 



[qn,Pn] = i^, on the conjugated variables of the eigen 
modes of the closed cavity. These commutation relations 
obviously translate to the commutation relations for the 
resonant variables, [QmPn] ~ i^, because of the canoni- 
cal nature of the transformations made in Sec. lIII Cl The 
unitary operator, which explicitly defines the correspond- 
ing quantum canonical transformation is 



U{t) = exp (^-iE(9' +Pl)^t/^fij 



(75) 



Averaging over rapid oscillations leads to the quantum 
Hamiltonian coinciding with the one in Eq. (24) with 



quantum operators replacing respective classical vari- 
ables. 

The quantization of the fundamental mode oscillator 
implies the quantization of the variable A{t) in terms of 
the conventional commutation relation for the annihila- 
tion operator, [A{t),A^{t)] = 1. 

Due to the linear coupling of the cavity to the trans- 



mission line, Eqs. 
ism outlined in Sec 
quantum regime 



the input-output formal- 



TV^ straightforwardly extends to the 
To this end, the classical amplitudes 
of the transmission line modes are to be replaced with 
the bosonic annihilation and creation operators, with 
[Ak{to) , (to)] = S{k — k'). From these commutation 
relations follows the commutation relation for the incom- 
ing field operator, [-B(i), {t')] = 6{t — t'), and similarly 
for the outgoing field operator C{t). 

The scattering relation, Eq. (41 1, has the same form in 
the quantum regime. 



C{t) = B{t) - iV2r A{t) 



(76) 



while the quantum Langevin equation for the cavity op- 
erator A{t) becomes, 

iA + dA + eA^ + a(^U + 1)A + iTA = %/2f B(t) . (77) 



This quantum Langevin equation, together with Eq. ( 76 ) 
preserves the commutation relation [^(i), A^^(t)] = 1 for 
the cavity mode, as shown in AppendixjD] The conserva- 
tive part of Eq. ( 77 ) is a dynamical equation associated 
with the Hamiltonian, 



H=-hs(A^A+l)-'^{A^^ 



^{A^A- 



-A' 



(78) 



- hV2riBA'' + sU) , 



B. Small quantum fluctuations 



A. Quantum Langevin equation 

The Hamiltonian description of the cavity parametric 
dynamics is a convenient starting point for the extension 
to the quantum regime. To this end we revisit Eq. (171 of 



Sec. Ill B and impose canonical commutation relations. 



The full analytical solution to the nonlinear quantum 
equation Eq. ( 77 1 is unknown. In what follows we restrict 



to the limit of small quantum fluctuations around the 
classical stationary states. Such a restriction is valid far 
from the bifurcation points and the parametric thresh- 
old. Some exact results for the critical fluctuations at 
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such points can be found in literature^SSSl^ also quantum 
jumps in multistable regimes have been investigatecP'^. 

To study quantum fluctuations within the framework 
of a Unearized quantum Langevin equation, we assume 
the cavity field operators to be of the form, A{t) — 
Aq + A{t), where Aq is a steady state solution of the 
classical nonlinear equation, Eq. (39), and A describes 
small quantum fluctuations, 



(79) 



Similarly, we separate the classical amplitude and quan- 
tum fluctuations of the input fleld in the transmission 
line, B{t) = Bo{t) + B{t), {B{t)) = 0. Then we expand 
Eq. ( 77 ) around Aq up to linear order in the quantum 



fluctuation A to obtain, 

ii + + ei^ H 
C = C + a\Ao\^. 



iVA 



2TB, 
aAl 



(80) 



Herein we introduced the effective detuning C, and the 
(complex) pump strength e by adding the terms propor- 
tional to the classical amplitude Aq. We note that Aq 
itself depends on the bare parameters 5 and e. Quanti- 
tatively, the parameter regions where this approximation 
is valid are identifled in Appendix [E] 

The analysis of Eq. ( pO| goes along the lines of 
Sec. |VD[ where the response to a classical detuned signal 
was evaluated. By introducing Fourier harmonics of the 
quantum fluctuations in the transmission line. 



B(4) 



dt 



B(t)e 



i<5fct 



and similarly in the cavity, the solution of the linear 



Eq. ( 80 ) is cast into the form. 



where 



A(<5fc) 
A\~6k) 



A- 



2VA- 



B{Sk) 
BH-Sk) 



C - 4 - ir -i 
-e* C + 4 + ir 



^(4) =(r-i4)' + C' 



l^(4)P = ( 



p2 ^^2 



Si 



(82) 

,(83) 
(84) 



2x2 



k ' 



cf. Eq. (60). It follows from this equation, that modes 
with frequencies Sk and —Sk are coupled pairwise by 
virtue of the parametric pumping. This property under- 
lines the generation of correlated pairs of photons with 
frequencies wi+a;2 = which is analogous to the photon 
generation under non-degenerate parametric resonance. 
The denominator in Eq. (82) turns to zero at 4 = 0, 



if the relation -f ~ I^P = holds, leading to the di- 
vergence of fluctuations at the corresponding parameter 
values. This happens at the parametric threshold, and 
at the bifurcation points, and indicates the enhancement 
of critical fluctuations. 



C. Fluctuations in the cavity 

The full power spectrum of the fleld in the cavity con- 
sists of the sharp line of the amplified (or generated) clas- 
sical signal, 27r|Aop^(4), together with the noise power 
spectrum, na{Sk), 



naiSk) 



dS',{AHS,)A{S',)) 



(85) 



Solving Eq. ( |82[ ) and assuming thermal noise in the input 
fleld, (^BH6k)B{S'^)) - N{Sk)5{Sk - S'^), where 7V(4) = 
^^h{n/2+Sk)/kBT _ 2) we calculate for the noise power 
spectrum 



naiSk) 



2F 



{|ef [iV(-4) + l] 



\D{SkW 

+ (F2 + (C-4)2)iv(4)} 



(86) 



At zero temperature, the noise power spectrum reduces 



to 



nnSk) 



2F|? 



\D{SkW 



(87) 



which can be interpreted as the amplifled vacuum noise of 
the input, manifesting itself as real photons in the cavity. 
The noise power spectrum in Eq. ( 87 ) has a resonance 



structure equivalent to the resonances in the classical re- 
sponse to a detuned signal discussed in Sec. |VD[ The 
only difference is that now the effective pump parameters 



enter Eqs. (83)-(84) instead of the bare pump parame- 



ters, since we allow here for a finite classical amplitude 
Aq. Accordingly, a single resonance a.t 6k — is observed 
under the condition 



< F^ + 



and otherwise two resonances are found at 



4 = ±VC 



F2 



(89) 



cf. Eqs. (63) and (64 1 



In Fig. 13 the noise power spectrum n.™'^(4) is pre- 
sented as a function of the pump detuning S for e = 2F 
and Bq — 0. In the monostable regime, S > — T^j 
where Aq = 0, the effective pump parameters in Eq. (80) 



are identical to the bare parameters, while in the bistable 
regime, S < \J — F^, with given by Eq. (45), they 
are C = 2Ve2 - F2 - 8 and jep = 5^ + r^. The condi- 



tion ([88]) identifies the interval {<? - 3FV2) js/t^ - < 



S < + around the parametric threshold, where 
the resonance lies at Sk = 0. Outside that interval, 
once the resonance is split, the separation grows with 
the parameter distance from the threshold, both below 
and above the threshold. At the parametric threshold 
itself, 6 = y/ e'^ — F^, the noise power diverges. 
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FIG. 13 . Noise power spectrum n™((5fc)-r of the cavity field, 
Eq. (871, vs. pump detuning 5. For 5 > \/(^ — V^ the classical 
cavity amplitude is Ao = 0, while |Aop > for & < — V^ 
according to Eq. ( |45[ ). The resonances, Eq. ( |89[ l, are indicated 
by the grey lines. (e = 2r, Bo = 0, a = T/lOO). 



The total number of photons in the cavity at zero tem- 



perature is {A^{t)A{t)) = l^op 
power 



with the noise 



d4 
27r 



l;|2 



/2 



r2 + c2 



(90) 



This quantity enters the validity criterium for the lin- 
earized Langevin equation, Eq. (79), which is analyzed 
in Appendix [E} 



D. Fluctuations of the output field 

Similar to the in-cavity field, the full power spectrum of 
the output field consists of the sharp line of the amplified 
(generated) classical signal, 27r|Cop(5((5A;), and the noise 
power spectrum ndSk), 



n-dSk) 



dS'k{cHSk)cis',)) 



(91) 



The relation between the input and output field opera- 
tors is similar to the one for a detuned classical signal in 
Sec.lVDl Eqs. (pl)-([59l). 



/ C(4) V ^ B{6k) 



(92) 



with matrix elements now dependent on the effective 
pump parameters. 



wi2 = 2irg, (93) 



and D{Sk) is given by Eq. (84). The matrix elements 
obey the fundamental relation. 



\iii{Sk)\^ - \vi2iSk)\^ = \DiSk)\\ 



(94) 



which provides the correct commutation relation for the 
output operators, [C{Sk),C^Sl)] = S{Sk - S'^). 



Equation ( 92 ) describes an input-output relation for 
a linear non-degenerate amplifiei*^ with signal and idler 
modes having frequencies Sk and — (5^, respectively, while 
the input classical tone at 6k = plays the role of an addi- 
tional pump. Indeed, the renormalization of the generic 
pump parameters in Eq. ( 80 ) is an effect of this addi- 



tional pump that increases the overall pump strength 
by cx a|^oP, and also affects the detuning 6 similar to 
Eq. ((28l). We note that Eq. (|92| is valid both below and 



above the threshold, and in the latter case it includes 
the classical parametric radiation acting as an additional 
pump signal even in the absence of the classical input. 

With the corresponding renormalization of the quan- 
tity q{Sk) = 2er/D{5k) that characterizes the amplifier 



gain, Eq. ( 62 ) , we cast the input-output relation, Eq. ( 92 ) 



into the form. 



(7(4)= + l9(4)Pe"^(*'='S(4) + ig(4)St(-4) 
= e'''('^'=) (coshri3(4) +sinhre'Xi?t(_4)^ ^ (95) 

where we introduced the standard notation for a non- 
degenerate parametric amplifier, 

sinhr(4) = |g(4)|, x(4) = argq(4)-?7(4) + - . (96) 

The mapping in Eq. ( |95[) is p rovided by a unitary two- 
mode squeezing operato^^'^, 

C(4) = e"^(''=)5K]B(4)5^K], (97) 



where ^ 



This implies that the stationary state of 



the output field is a pure state provided the input is a 
pure state. This is true in spite of because the evolution 
of the total system, including the cavity variable, is for- 
mally non-unitary due to the presence of the dissipative 



term in the Langevin equation ( 80 ) 



The noise power spectrum of the output field can be 
computed from Eq. (95), and for thermal noise input it 
reads, 

n,(4) = N{5k) + |g(4)|' [N{5k) + Ni-5k) + 1] .(98) 
At zero temperature this equation reduces to 

nr(4) = |g(4)|' = 2rnr(4), (99) 

and describes the generation of real photons from the 
vacuum under parametric resonance. This phenomenon 
is closely related to the Dynamical Casimir effect - the 
creation of real pho tons from vacuum fluctuations by 
an accelerated mirroi!2Sl22l. Here the role of the moving 
mirror is played by the time-dependent boundary condi- 
tion, driven by the modulated magnetic fiux through the 
SQUID. 

The output noise, being proportional to na{6k), inher- 
its the resonant behavior of the noise power spectrum 



in the cavity, as discussed in Sec. VII C and shown on 
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le+06 



, 2 
S/T 



6e-05 



FIG. 14. Noise photon flux n^Vr, Eq. fToTj l, vs. 5 and e, 
assuming (a) bare pump detuning, and (b) taking into account 
the pump-induced frequency shift, Eq. (281. (Bo = 0, q = 
r/100). 



Fig. [TSj In the deep subthreshold regime, for very weak 
pump strength, e ^ F, and in absence of an input signal, 
Bo = 0, Eq. (Iggl) takes the form. 



=(4 



4e^F' 



2-r2 



[F2 + (4 + Sr] [F2 + (4 - 6)^] 



(100) 



In this limit the resonances move tovifards Sk — ±(5, and 
the resonant structure of n™'^(4)i approaches the one 
computed irP^l and observed irf23. 

The total photon flux in the output field is 



{CHt)Cit))^\Co\' + nl 



with the noise photon flux, 
Fle-P 



F2+C2-|e|2 



(101) 



at zero temperature. Below the parametric threshold, 
e < y/S^~+T^, the effective parameters in Eq. (80) 



are identical to the bare ones, and Eq. (101) reduces 
to n™'^ = Te^/ {T^ + - e^). Above the threshold. 



e > V4 +F2, with |Ao| > given by Eq. (|45]), Eq. ( [lOl| ) 
becomes 



F((52 + r2) 



4Ve2 -F2(-5 + 



F2) 



(102) 



The output noise level is illustrated in Fig. [14] as a func- 
tion of e and S for Bq — 0. The right panel demon- 
strates the effect of back-bending of the threshold line 
due to the pump- induced frequency shift, Eq. ( [28| (cf. 
Fig. [gb) in Sec. |VA|. The noise is enhanced at the 
parametric threshold and decreases while moving away 
from the threshold, there it is estimated as n™'^ ^ F for 
e, (5- F. 

Since the noise near the parametric threshold becomes 
strong, it is useful to evaluate the conditions for the out- 
put coherent signal dominating over the noise, |CoP >• 
(CtC). 

Above the parametric threshold, the signal-to-noise ra- 
tios are identical for the output field and the field inside 
the cavity (for Bq = 0), 



|Co| 



l^ol 



(103) 



Therefore the limitation established by Eq. (E3) for the 



field in the cavity applies as well to the output field 

F /a\2/3 



F> 



(104) 



Below the threshold, the maximum amplified signal is, 
according to Eq. (56), 



|Co| 



1 + + q \Bo\' = 



(r- 



(F-6)^ 



\Bn 



(105) 



for ^ = and a|AoP ^ \/F2 — e'^. Comparing this with 



Eq. ( 101 ), we arrive at the constraint on the input signal. 



\Bn 



> 



e2 1-e/F 
F (l + e/F)3 



(106) 



This bound is of order F for e ^ F, and decreases both 
at weak pumping and close to the threshold. This is 
explained, at small e ^ F, by the fact that the amplifi- 
cation of vacuum noise is small, while the classical signal 
remains finite, and, close to the threshold, by the fact 
that amplification of the signal is more efficient than the 
amplification of the noise. The constraint in Eq. ( 106 ) 



is qualitatively similar to the one for the field inside the 



cavity given by Eq. ( E6 ) 



E. Squeezing 

A homodyne detection scheme allows for measurement 
of the quadratures of the outp ut signa l, and characteriza- 
tion of quadrature fluctuationJ3IIllIlll_ With this method, 
the output field is mixed with a strong classical field of a 
local oscillator, B^o cos{flt/2 — 6), and the intensity of 
the mixed signal is measured. This intensity is propor- 
tional to the output field quadrature, loit) ~ BlqX^ {t): 



(107) 



The phase 6 refers to the phase shift of the local oscilla- 
tor with respect to the parametric pump; variation of 9 
allows accessing all the quadratures individually. 

The mean quadrature is determined by the classical 
output signal 

{X'^) = = Coe-"' + C*e"' = 2|Co| cos(0c-^^) • (108) 

Separating the classical and quantum components, 
X%t) = X^+X\t), X\t) = (7(t)e-"'+C't(f)e'^ and us- 
ing the spectral representation of the noise quadratures, 
X%6k) = C{Sk)e-''^ + C^-Sk)e''', we present the cor- 
responding power spectrum in the form 27r(Xg )2|5(i5fc) + 
5^(4), where 



^^4 



dSUx<>{6k)xOiS',) 



(109) 



is the squeezing power spectrunP^EH. Note that by virtue 
of the stationary state of the cavity, (X^{6k)X^{S'i.)) oc 
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6{5k+5'^.), hence only symmetric correlations between the 
sidebands contribute to the integral, i.e. the squeezing 
power characterizes the two-mode squeezing. 

We calculate the squeezing power assuming vacuum 
fluctuations of the input, using Eq. (95) for the output 
field operators. The result reads. 



2\q\' - 2v/rTi^Im [q*{-Sk)e'^^~^'^ 

4r 



Im(?e-2i9 



2r|e|2 

r2_ 



)( 



I- 2rCRe(ee"2'^ 



Si 



(110) 



Equation S^{Sk) = 1 corresponds to pure vacuum fluctu- 
ations. The noise squeezing power varies with the phase 
0, the maximum and minimum values reached at and 
^0 + T^/"^! respectively, with 



tan(26'o) = 



2Crim(?) 



[T^-e + \i\' + dl]ReCe) 



2CrRe(e) + [r2 - + |e|2 + 52]Im(e) 



The corresponding extreme values are determined by the 
quantity |(7((5/i;)|, and have the form. 



controlled by the input signal phase 6b- By varying the 
latter one may control the signal-to-noise ratio for the 
quadratures 
In Fig 

0, and two different values of e > F 



for Bo = 



15 'a-b) the squeezing power S (5^) is shown 



The 



^-values of maximum and minimum squeezing power are 
indicated by black lines. 

It is useful to also quantify the in-cavity squeezing, 
by calculating the squeezing power for the quadrature 
operator (4) = A(Sk)e~'^ + A\-6k)e'^ in analogy to 



Eq. ( 109 ). Although the phase 9 in this case is not related 



to any externally tunable phase, it might be relevant for 
a quadrature-dependent coupling to a qubit placed in the 
cavity, or to another transmission line. 

Assuming vacuum fluctuations in the input field we 



calculate the internal squeezing power using Eq. ( 82 ) 



2F 



+ (C - Skf + kf - 2CRe(ee-2'^) 

(113) 



\D{6k 

-2(F- i4)Im(ee"2' 



Further evaluation of the minimum uncertainty of the 
cavity quadrature, ((AXf)^) = (l/27r) / (4), re- 

sults in the value 1/2, as in the case of linear parametric 
amplifierJ^, i.e. a factor 1/2 below the vacuum limit. 
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FIG. 15. Squeezing power of output noise 5* Eq. (110 1 

(a-b), and in-cavity noise |S'®(5fc)| ■ F, Eq. (113 1, (c-d) for 
e = I.IF (left column) and e = SF (right column). The 
quadrature phases for maximum and minimum squeezing 
power, Eq. ( |111| are indicated by the black lines. {5 — 0, 
Bo = 0, a =¥7100). 



2 ' 

5«-«''+-/2(4)^ (yrTM4yF±i9(4)i) ,(112) 5^' 

which is similar to the classical gain of the ideal amplifier, 
Eq. ([57]), including the relation, 5^o(4)-5'^''+''/^(4) = 1- 
However, the maximum squeezing and maximum quadra- 
ture gain do not generally correspond to the same value 
of mixing phase 9. Moreover, the amplified classical sig- 
nal, Eq. (1081, contains an additional phase, Oc, which is 



igs. [Tslc-d) 15*^(4)1 is illustrated for = 0, (5 = 
e > 1 , in comparison to the external squeezing 



In Fi, 
and 

5« (4) of Figs 

detuning in Eq. (113), >S'^(4) is not symmetric around 
4 = 0, as is the case for 5^(4)- 



15 



a-b). As a consequence of the effective 



F. Second order coherence 

The two-mode squeezing is a nonclassical property of 
the amplified noise that originates from the production 
of noise photons in entangled pairs. Further information 
about the nonclassical properties of the correlated output 
photons is provided by a two-photon correlation function, 
and characteristics of two-photon entanglement. 

We start with evaluating the second-order correlation 
functiorPsl, 



G(2)(r) = {C\t)C\t + T)C{t + T)C{t)) . (114) 

In the presence of the classical output component, this 
equation takes the form, 

G(2)(r) = \Co\^ + 2\Co\'{{C\t)C{t + T)) + {C\t)C{t)) 
+Re(e-2iec(c'(i)c'(i + ^))^| 

+ {C\t)C\t + T)C{t + T)C{t)) . (115) 

Explicitly, using Eq. ( |95[ ), we obtain for t — Q and input 
vacuum noise, 



G(2)(0) = (|CoP 



(116) 
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FIG. 16. (a) Normalized second order correlation ^'^'(O), 
Eq. ([TTel, vs. e for fixed Ob = -7r/4 and for S/T = 0,1,2 



(from bottom to top, the corresponding parameter values are 
marked with crosses on panel (b)). (b) g'^'(O) and (c) eam 
IColVl-BoP vs. S and Ob, for fixed e/V = 0.2 (indicated by 
the vertical line in (a)). (|Bo|^ = T/IO, a = T/lOO). 
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In Fig. |16[a) the normalized correlation function, 
(2)(0)^= G(2)(T)/(Ct(i)CW>' = G(2)(r)/(|Co|2 + 
n™'^)^, is presented as a function of the pumping strength 
e for several values of the pump detuning S. In the Duff- 
ing limit, e = 0, all the terms in Eq. (1161 vanish ex- 



cept of the first one, yielding the coherent state limit, 
(7^^-'(0) = 1. The same is also true for large pumping 
strength above the threshold, e > + 5^ This is ex- 
plained by the rapid growth of classical radiation power 
that dominates over the fluctuations, |CoP ^ n™*^ (kinks 
on the curves at e/F > 1). 

At the intermediate pump strengths both bunching 
(g(2)(o) > 1) and antibunching (g'-^^O) < 1) are pos- 
sible. For pure output noise in the absence of classical 
output, Co ~ (i.e. for Bq = below the threshold), 
only bunching occurs, ^^^•'(O) = 2 + {T'^ -I- C^)/|ep, where 
the degree of bunching exceeds that of classical chaotic 
radiation, (7(^^(0) > 2. This can be interpreted as a con- 
sequence of the pair production of noise photons. 

When Bq > 0, also antibunching is possibl^^ due to 
the interplay between the classical and the quantum con- 
tribution to the correlation, last line in Eq. (116). It 



occurs within a relatively narrow window of parameters, 
e < F, (5 > F, |i?oP ~ F, for which the phase dependence 
in the last term in Eq. ( 116 1 can introduce a sign change. 

The dependence of .g^^^(O) as a function of the in- 
put phase 6b and the pump detuning S is illustrated 
in Fig. |16[ b). Pronounced antibunching (blue regions) 
is observed for |i3oP/F < 1, and for those values of 6'^ 
where the gain approaches unity, |CoP/|i?oP ^ 1, com- 
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FIG. 17. Entanglement entropy E (solid) and logarith- 
mic negativity Ej^/ (dashed) vs. pump strength e for 5k /T = 
0.02,0.2,1.0,2.0 (from top to bottom). The vertical lines 
mark the values of e at which the resonance, Eq. (891, is 
encountered at the chosen value of Sk- {S = 0, Bq — 0, 
a = r/100). 



pare Fig. |16[c). 



G. Two-mode entanglement 

The degree of entanglement between the two modes 
with frequencies Sk and can be quantified with the 
entanglement entropj0^, 



i;(4)--Tr(p(Jfe)lnp(4)) 



(117) 



where p{Sk) is the reduced density matrix of one of the in- 
volved modes. If these modes are entangled, the entropy 
takes a positive value, E > 0. 

We compute the entanglement entropy for the ampli- 
fied vacuum noise, using the two-photon wave function 
of the squeezed state, 

00 

\6k,-Sk) = (coshr)-i^(tanhr)"e™x|n,n) , (118) 

ji=0 

which is obtained by applying the squeezing operator, 
Eq. (97), to the vacuum input, \6kj—Sk) = S[^] |0), and 
using the decomposition equatiori-9'. The reduced density 
matrix has the form, 

CO 

p{5k) = (coshr)-2 ^(tanhr)2"|n)(n| , (119) 

11=0 

giving the entanglement entropjBH^ 

E = cosh^ r ln(cosh^ r) - sinh^ r ln(sinh^ r) . (120) 

The entropy is nonzero for all e > 0, and follows closely 
the squeezing parameter r{Sk), asymptotically approach- 
ing the linear dependence, E ^ 2r, for r > 1. 

The entanglement entropy E is shown as function of e 
in Fig. 17 (solid lines), for several values of the detuning 
6k, and for Bq = and (5 = 0. For small detuning, 5k < 
F, the entropy reaches the maximum at the threshold. 
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e = r, at which r{5k = 0) = 0)] diverges. With 

increasing value of the detuning Sk this maximum shifts 
towards the value of e((5/j), at which q{6k) exhibits the 
resonance, Eq. (89). The entropy rapidly decreases above 
the threshold, analogous to the behaviour of (7(^^(0), due 
to the emergence of the classical radiative state, Aq ^ 0, 
that suppresses q{Sk)- 

A convenient measure of entanglement for Gaus- 
sian states is provided by the logarithmic negativitj^^ 
related to the covariance matrix for the two entan- 
gled modes. The covariance matrix is defined 
through a 4-vector composed of the quadratures, — 

14/3 = I {RaRp + RpRa) - {Ra) {Rfi) ■ (121) 

Then splitting the covariance matrix into 2x2 subma- 
trices, V = (V^i, V3; , V2) , the logarithmic negativity is 
defined as 



Ej^ = max (0, — ln(;/_)) 



(122) 



where = y (cr - Vcr^ - 4 det V) /2 and cr = det Vi + 
det V2 — 2detV3. For entangled states the logarithmic 
negativity takes positive values. 

For amplified vacuum noise we obtain a simple result, 
using Eq. (95), 



E^f = 2 In 



yiTMW + Wk)\) = 2r , (123) 



i.e. the logarithmic negativity is equal to twice the 
squeezing parameter r{5k). 

The logarithmic negativity Ej^ is shown in Fig. |17| with 
dashed lines. Its functional behavior is basically equiva- 
lent to that of the entropy E. 

Our calculation shows that the degree of the two-mode 
entanglement is significantly enhanced in the presence of 
the parametric resonance. To evaluate the exact max- 
imum entanglement value one needs to go beyond the 
quasilinear approximation and include the nonlinear ef- 
fect. We make a qualitative estimate by taking the func- 
tion |(?((5fc)l the threshold, e = T , 8 — and at 5k — 0, 
and for the cavity field given by Eqs. (66) and (68) as- 



suming input power, I-BqP ~ F, corresponding to one 
photon per bandwidth. This yields an estimate. 



Er, 



const • In ■ 



(124) 



with a numerical constant of order one. This crude es- 
timate seems to agree with more accurate evaluation of 
the critical fluctuations^S. For values F/a ~ 100 achiev- 
able in tunable cavities, the entanglement entropy can 
accordingly reach the values 4.5 — 5. This is significantly 
larger than the values calculated^ for a non-resonant 
open transmission line with modulated boundary, and 
also exceeds the values repor ted for experimental para- 
metric Josephson dey'ice^^^. 



VIII. SUMMARY 

We have developed a consistent theory of parametric 
resonance in a high quality tunable superconducting cav- 
ity. We considered the nonlinear classical dynamics of the 
cavity both below and above the parametric threshold, 
and analyzed amplification of external signals, and para- 
metric radiation. We also studied quantum properties of 
the amplified and radiative fields. 

The non-equidistance of the cavity frequency spectrum 
enabled us to formulate the theory of the degenerate 
parametric resonance in terms of the one encountered 
in a nonlinear parametric oscillator. We identified the 
parameters of this effective oscillator as functions of the 
cavity generic characteristics, and investigated the multi- 
stable cavity dynamics in a relevant range of the effective 
parameters. 

The operation of the device in the monostable regime 
as a nonlinear parametric amplifier is characterized with 
a phase-dependent differential gain, which increases at 
small input power and reaches the maximum value at 
the parametric threshold. We found that this maximum 
value scales with the ratio of the damping coefficient and 
the nonlinearity coefficient, F/a. We also found that the 
relation between the maximum and minimum gain for an 
ideal linear amplifier is violated in the nonlinear regime, 
Gmin ^ 1/Gniax- Extremely small values of a available 
in tunable cavities allows for very large gain and strong 
amplification vs. deamplification contrast. 

Amplification of detuned signals was found to exhibit 
sideband resonances within a specific region of the cavity 
parameters. This effect can be used for enhancing the 
amplification bandwidth while maintaining high gain. 

The application of the device as a parametric bifurca- 
tion amplifier was discussed in regard to dispersive qubit 
readout. The advantage of the parametric regime com- 
pared to the conventional JBA is a high sensitivity of 
the strength of the output signal to the variation of the 
cavity frequency. This, together with a high amplifica- 
tion gain, provides a potential for improving the fidelity 
of qubit single shot readout. 

Yet another suggested method for qubit readout is 
based on a high contrast between the strengths of para- 
metric radiation above the threshold and amplified noise 
below the threshold. 

Small-amplitude quantum fluctuations around the 
classical signal were investigated for the in-cavity field 
and the output field. The limit of small fluctuations is 
appropriate in a wide range of the device parameters ex- 
cept of small regions of critically enhanced fluctuations 
close to the bifurcation points and the parametric thresh- 
old. The theory is analogous to the one for a quantum 
linear amplifier. The strength of the amplified noise in- 
creases in the vicinity of the threshold in accord with the 
classical gain. The same is also true for the two-mode 
squeezing and the entanglement quantified with the en- 
tanglement entropy and the logarithmic negativity. At 
the threshold, the estimated magnitude of the squeezing 
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parameter may reach the values of a few units, exceeding 
that achievable e.g. in non-resonant Josephson mixers. 

The second order coherence is dominated by strong 
bunching for small classical inputs, resulting from the 
production of noise photons in pairs. However, for classi- 
cal inputs with strength comparable to the vacuum noise, 
significant antibunching is predicted resulting from the 
interference of the classical and quantum field compo- 
nents. 

To conclude, we note that the developed the- 
ory straightforwardly extends to the regime of non- 
degenerate parametric resonance, when the pumping fre- 
quency is commensurate with a combination of cavity 
resonances. Similarly, in this case, strongly enhanced 
amplification gain is to occur near the parametric thresh- 
old, as well as strongly enhanced two-mode squeezing and 
entanglement of the cavity modes selected by the reso- 
nance. 

Yet another extension of the theory is readily done 
for a two-sided cavity parametrically pumped by two 
SQUIDs, attached to both sides of the cavit^. The 
dynamics of this device is equivalent to the single-sided 
parametric cavity, provided the SQUIDs are operated 
at the same pump frequency. The parametric resonance 
is then controlled by an effective pump strength, which 
depends on the phase shift between the actual pumps. 
For equal pump amplitudes the parametric effect is max- 
imum for the out-of-phase pumping ( "breathing" mode) , 
while for the in-phase pumping ( "translational" mode) 
the parametric instability is completely suppressed. 
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Appendix A: Lagrangian and boundary condition of 
the flux-tunable cavity 

In this appendix we derive the Lagrangian of the flux- 
tunable cavity, Eq. ([2]), and give arguments for its valid- 
ity. 

We start with a description of the SQUID establishing 
the connection between the cavity and the pump line. 
The generalized coordinates of the SQUID are the su- 
perconducting phase (j)d = (j^i^ = d,) at the cavity edge 
X = d, the phase 2/ dropping over the inductance L of 
the SQUID loop, and the phase 2/oxt dropping over the 
coupling inductance L^xt of the pump line, see Fig. [T] 

The SQUID is modelled as symmetric, with two iden- 
tical Josephson junctions, each having a Josephson en- 
ergy Ej and a capacitance Cj. To simplify notation 
we assume that the SQUID is grounded in such a way 
that its geometric inductance L is divided into two equal 
parts L/2, with a phase drop of / over each part. Thus, 
the phase difference on one of the Josephson junctions is 
0d — / and (t>d + f on the other. The coupling to the flux 



line is inductive, with a mutual inductance M <^ L, L^xt- 
The full SQUID Lagrangian is 

2 

-/) + i;jcos(0d + /) 



h 

-Ej cos 

h 

2e 



/ + 



2 



1 

2L 



4/^ 



8 //oxt 
J^ext 



(Al) 



or, written with the capacitive energy of a Josephson 
junction Ec = (2e)^/(2Cj) and the inductive energy of 
the SQUID loop El = {h/2ef/{2L) 

' ^ - 2Ej cos {(pd) cos / 



2Ec 



2E, 



c 



4/2 



M 



-ffc: 



(A2) 



Separating the (/)d-dependent terms (first line) from the 
purely /-dependent ones (second line), £s = £s[4'd, f] + 
Cs[f], the former can be combined with the bare cavity 
Lagrangian 



£(0) 



dEr 



2^2 



dx 



(02 2^,2) 



(A3) 



with the inductive energy of the cavity -EL.cav — 
/ {Lfjd{2eY) . Together, these form the Lagrangian £cav 
of the flux-tunable cavity, Eq. ([2]). 

For typical cavity and junction dimensions the or- 
ders of the three inductive energies in the Lagrangian, 
Eqs. (A2)-(A3), are distinctly different. The dominant 
energy, Eijh^ IQ^ GHz, determined by the small ge- 
ometric inductance of the SQUID loop (L w 10^^2h), 
is larger than the Josephson energy of the SQUID, 
2Ej/h ^ 10'' GHz, and that dominates over the in- 
ductive energy of the cavity, E^^cav/fi ^ 400 GHz (for 
dLo f« 2 • 10~^H). Furthermore, the Josephson plasma 
frequency ujj = \/2EjEc /fi. ~ 300 GHz is high compared 
to the fundamental cavity resonance, ujq ^ 40 GHz (com- 
pare Sec. HI A). 

The equations of motion for (j)d and /, according to the 
full Lagrangian £cav + ^sif], 

■■ 

— (f>d + '^Ej cos/ sin -I- EL.ca.YC 
Ec 

■■ ( 

-f+Ej cos ^dSinf+AEL / + 



2Er 



i>d^O 

M/cxt 



(A4) 



= (A5) 



describe two coupled nonlinear oscillators. For /oxt — 
the equilibrium is (/ = 0, = 0). 

In general, the coupled dynamics of nonlinear, driven 
oscillators features chaotic behaviour. We restrict our 
analysis to the case (fid ^ and assume that this is ful- 
filled even in the presence of a resonant excitation by 
the external field /oxt(0- Under this condition the equa- 



tion of motion for /, Eq. (A5), decouples from the other 
oscillator, 

fi2 .. f M \ 

/ + Sjsin/ + 4Si /+— /ext =0, (A6) 



2Er 
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and the dynamical equation for (pd, Eq. (A4), then de- 



pends only parametrically on /(i), of. Eq. ([7|) 

We suppose the external force of the form /oxt(0 = 
Fext + 3fcxt{t), Sf cxt ^ 1, and separate the SQUID phase 
response f{t) — F + Sf{t) into a constant equilibrium 
shift F, governed by the equation, 

Ej sin F + 4El {F + {M / L,,t)Fe,t) = , (A7) 
and a small harmonic oscillation, 5f{t) <^ 1, driven by 

<5/ext(i), 

^Sf + (E J COS F + 4EL)Sf = -i^5/,,t(t).(A8) 

Assuming (5/oxt(i) = ^/cxtCosi7t, we write the stationary 
solution in the form, Sf{t) = Sf cos fit, 



SELMEc/h^L,,t „ 

A/ = 2 7^ '^/oxt , 



(A9) 



where ujj ^ uj \J cos F + A,El /Ej is the frequency of the 
/-oscillator, which is much larger than the frequency of 
the pump, VL « 2^0- 



Linearized around the equilibrium shift F, Eq. (A4) 
becomes 



Er 



+ 2Ej [cos F — sin F5f cos{Q.t)] sin (jid 
+ EL,c^^d(l)'d = . 



(AlO) 



For 6f = this boundary condition determines the cavity 
mode spectrum, Eqs. ^ and ([9|. 

Further expanding Eq. (A4) to the second order with 



respect to Sf leads to the pump induced shift of the cavity 
frequencies. Indeed, averaging over time, we get a cor- 
rection to the Josephson energy, 2EjcosF{l — Sp/A), 
which will modify Eq. ^ accordingly. In particular, for 
the fundamental mode we get from Eq. (10), 



LOQ (Sf) 

Wo (0) 



(All) 



One could also expand Eq. ( A5 1 to the second order with 



respect to 6f, which would lead, after the time averag- 
ing, to a shift of the static bias F, and eventually to an 
additional shift of the cavity frequencies. However, this 
effect is small, by virtue of the parameter Ej/El ^ 1, 
compared to the shift (All). 



Firstly, making use of Eq. ^ , the overlap integrals of 
the non-orthogonal modes are 



dx COS kjf^x cos kjYi^x 



dM„, 



2Cj 



cos kn d cos km d 



dklM,„ 



dixkn ktn, sm k^i X sm k^yj X — ^ 



(Bl) 



(B2) 



2Cjkl 
Co 



cos knd cos kmd — km cos fc„d sin km,d , 



where we have defined the coefficients Mn, Eq. (14 1 



With these, the bulk contribution to the cavity La- 
grangian becomes 



2eJ 2 7„^K^^-«W 



(B3) 



2^1 Co 



cos knd COS kmdq 



2 / 2 \ 

+v COS knd ( -^fcm COS kmd + km sin kmd j gnQn 

In the remaining boundary contribution of Eq. ([2|, 
we firstly separate a time dependent, nonlinear poten- 
tial term 

V{cj)d, t) = -2Ej (^cos fit) cos + cosF^'^ (B4) 



from the harmonic contribution, 



h \ 2C 

^ -^^l + 2Ej cos fit) cos<j)d = 



(B5) 



The mode-representation, Eq. (12), of the harmonic part 
becomes, using i2eYlif?Co) = tr/(rf£'L,cav), 



hy 2Cj -n 05 
2Cj 



(B6) 



E cos knd COS kmd qnqm 

COS knd COS kmd qnq-n 



2Co 
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2Ejv^cosF 



Appendix B: Mode representation of cavity 
Lagrangian 

In this appendix, we express the Lagrangian of the 
fiux-tunable cavity, Eq. ([2|, in the mode representation, 
Eq. (13), based on the expansion (12) of the cavity field. 



The first term of this cancels directly with a term in the 
bulk contribution, Eq. ( B3 ) . Further, using the definition 
of the modes in Eq. ([9|, we note that 

fc, I + sin k..A = (B7) 
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leading to further cancellation of terms between the bulk 
and the boundary contribution. The remaining terms are 

-Ccav =IJ2 [^'^»^" - MnV^f'lll] - ^(ln,t) , (B8) 



with V{qn,t) = V{(l)d{qn),t). This is the mode represen- 
tation of the cavity Lagrangian in Eq. (13 1. 



line is evaluated (for a; < 0) 
dk2cos{ijJk{t — i')) cos{kx) 

dk [cos(fcu(t — t') — kx) + cos{kv{t — t') + kx)] 



dk 



'sm{kv{t - t') - kx) + sin(fcw(t - t') + kx)] 



Appendix C: Transmission line amplitudes and 
scattering relation 

In this appendix we show the relation of the flux am- 



dt Jq kv 

^ [sgn(t -t' - x/v) + sgn{t - t' + x/v)] 
Zv dt 

- [S{t ~t' - x/v) + 5{t -t' + x/v)] . 

V 



(C6) 



plitude B introduced in Eq. (40) to the incoming field, 



and similarly for the flux amplitude C of the outgoing 
field, as well as their mutual relation given in Eq. (41 1. 



The incoming and outgoing fields in the transmission line 
are defined, respectively, as 



'dk r 



, , , e nr f^dk r 

t>out{t)=^\ / ^= 

n V ttCo Jo y/ujk 



afe(io)e-'"'=(*-*''^+h.c.l (CI) 



The second i5-function gives a contribution at t' = t + 
x/v^ whereas the first, for x < 0, is not included in the 
integration limits, t' = t — x/v > t, and therefore 

dk / dt'pU') cos(ujk(t — t')) cosikx) ~ — y(t + x/v). 
10 Jto 2w 

(C7) 

Taken together, the field in the transmission line reads 

0TL {X,t) = (f)in{t - x/ v) + 4)in (t + x/v) 

n Cod V OqMo V 



afc(^i)e~''^'°'* +h c (C2) Alternatively, the transmission line field, Eq. (29), can 



be evaluated from the second solution for the ak(t), 
Eq. ([C4|, yielding 



These are based on the solutions of Eq. (34), which is 



expressed in terms of initial amplitudes akjto) at a time 
to < t in the past, or in terms of final amplitudes ak(ti) 
at a time ti > t in the future. 



(j)TL {X,t) ^ (f>out {t - x/v) + 0o„f {t + x/v) 

2e Cc I 2dujo 1 



h Cod V CqMo v 



p{t ~ x/v) . (C9) 



afc(i)=«fc(io)e-'"''(*-*°) 

Cc 1 2dujoujk /■* 
Cody MoTTh 4 

ak{t) ^ ak{h)e-'^''^*-'^^ 



dt'p{t')e 



-iu:k(t-t') 



(C3) 



(C4) 



By subtracting Eqs. (C8) and (C9| at x = we can es- 



tablish a relation between the incoming and the outgoing 
field components. 



Cc / 2dujouJk 



Cod\ Monh 



dt'p{t')e 



-iojkit—t') 



Such a definition is justified, as will be shown below, 
by attributing different propagation directions along the 
transmission line for the incoming and outgoing field 
components, which can be separated by circulators and 
hence have physical meaning. 



J. J. '2e Cc 2dujo p{t) 

n Cod V MoCo v 
Note that the last term can also be expressed by 
the derivate of the cavity field at a; = 0, (f>o{t) « 
{2e/h)y/2ujo/iCodMo)p{t), using p{t) w ujoq{t) in the 
weak coupling approximation. 



We can now evaluate Eq. (CIO) at t + x/v, and in 



serf in Eq. ([C8j), such that the transmission line field is 
expressed as a linear combination of (/)i„ and (f>out alone, 



{X,t) = (jiin (t-x/v) + (f>out {t + x/v) , 



(Gil) 



We firstly use the solution (C3) to evaluate the field in 
the transmission line, Eq. (29 1 with qk = \Jh/2(ak + 



demonstrating the role of and (/)out as incoming and 
outgoing field components. 

Finally, we want to relate the general input-output re- 



(l)TL{x,t) = 2e 



' dk cos kx r 

\/nhCoUJk 

t 



afc(io)e-'"''-('-*°)+h.c.l(C5) 



+ / dt'p{t') I dkcos{u}k{t~t'))coskx. 

TihCod V CqA/o Jto ^0 

The contribution from the first line can be straightfor- 
wardly identified with from Eq. (CI), and equals 



lation, Eq. (CIO), with the slow varying amplitudes of 
the resonant approximation. To that end we separate 
the fast time oscillation with frequency $7/2 in Eqs. (CI ) 
and (ICl, 



2e 

h V 2Cot<Jo 
2e 



- {Bit). 



Ant/2 



4'in{t — x/v) -\- 4>in{t -\- x/v). The integral in the second 



w(t) = ^/^^(c(t)e--/^-^h.c.),(C13) 



h.c. (C12) 
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where 
Bit) 
C{t) 



27r 



2tt 



dk 



V^fc 



afc(io)e~'"''(*-*«) (C14) 
afc(ii)e-'"''(*-*^).(C15) 



Within the resonant approximation, 



Wo, 



iOt/2 



Afc, these quantities coincide with the ones 



defined in Sec. IV cf. Eq. (40). The cavity momentum 



is expressed in the rotating frame as well, 



p{t) = -i 



A{t)e 



A\t) 



^iOt/2 



(C16) 



with the slowly time-dependent cavity amplitude A{t). 
By setting these expressions into Eq. (CIO I, multiply- 
ing with e'^*/^ and averaging over fast oscillation, the 
corresponding input-output relation is obtained in the 
rotating frame, cf. Eq. (41). 



C{t)-B{t) = V^ ^(^) ^ . (C17) 
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Appendix D: Commutation relations 

In this appendix we show that the quantum Langevin 
equation, Eq. (77), preserves the commutation relation 



of the cavity amplitude. To this end we express the so- 
lution of Eq. (77) in terms of the propagator U{t,to) — 
exp(— i_ffcav(^ — io)/^) with the Hamiltonian that governs 
the dynamics of the isolated cavity, 

i?cav = -M(^AU+i)-^((At)2 + A2) (Dl) 



ha 



(cf. Eq. ([78|)). The operator A refers to the Schrodinger 
picture and coincides initially with the Heisenberg oper- 
ator A{to) = A. At time t the solution is 

Ait) = e-r(*-*«)c/-i(t,to)A(to)C/(t,to) (D2) 
-iV2r f dt'e-^'^'-*'^U-\t,t')B{t')U{t,t') 

J to 

= e-^^'-*°^U-\t,to)A{to)U{t,t„) (D3) 
iV2r / dt'e-^^'-*">B{t'), 



where we have used the fact that Hav and B{t') com- 
mute since A is uncorrelated with the operators ak(to) 
of the incoming transmission line modes of which B{t) 
is composed. Using this solution we are able to evaluate 



the equal time commutator, 

[Ait),A^{t)]^e-^^^'-'«^[U"\t,t,)Aito)U{t,to)Am 
U-\t, to)AHto)U{t, to) ] + iA/2fe-r(*-*«) 



X dt'e 
It, 



2T 



i-Tt 



' {[U-\tMA{to)U{tM), B\t')\ 
- [B{t'), U-\tM)AHto)U{tM)] 
dt'dt"e~^^''+*"'^ [B{t'), B^t")] . (D5) 



The mixed commutators vanish, again with the ar- 
gument of initially uncorrelated cavity and transmission 
line operators, leaving 

[Ait),A\t)] = 

e-^^^'-'°^U-\t,to) [Aito), A\to)] U{tM) 

+2T [[ dt'dt"e-^^''+'"^ [B{t'), B\t")] . (D6) 



to 



Finally, using [A{to), A^t^)] = 1 and [B{t'), B\t")] = 
5{t' ~ t"), we arrive at the desired result, 

[Ait),AHt)] = e-2r(t-to)_(^e-2r(t-t„) _ i) = i . (D7) 
The invariance of the commutation relation under the 



Langevin evolution, Eq. (77 1, follows from the correct 



combination of the damping term, TA, and the fluctu- 
ations in the amplitude B. Averaging over fluctuations 
would violate the exact unitary evolution and break the 
commutation relation. 



Appendix E: Validity of quantum linearized 
treatment 

Having evaluated the magnitude of the quantum fluc- 
tuations, we are able to discuss the region of validity 
of the linearized equation, Eq. (80 1. Two assumptions 



have been made for the derivation: the amplified signal 
at frequency 6k = has been treated as a classical field, 
|Aop S> 1, and its magnitude to exceed the amplified ex- 
ternal noise, |AoP ^ rt^. Together these conditions are 
(cf. Eq. (l79|). 



\Aq\^ > max(l, Ua) 



(El) 



We analyze these conditions separately above and below 
the threshold, at zero temperature, and at 5 = for 
simplicity. 

Above the threshold, e > F, the parametric radiation 
dominates over the input signal. Neglecting the input, 
Bq = 0, we have \Ao\'^ ~ \J ~ /a in accord with 



Eq. (|45l. Then |e|2 = F^ and C = 2y/e^ -V^, and the 



amplified vacuum noise is. 
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The conditions of Eq. (El ), 
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r 2V2 



(E3) 



are fulfilled everywhere except of the close vicinity of 
the threshold. Near the threshold the external noise 
dominates, while its role diminishes with growing pump 
strength. In the limit of very strong pumping, e 3> F, 
Eq. (E3| reduces to a/e <C 1. This result can be under- 
stood from a purely Hamiltonian argument. The semi- 
classical limit requires the quantum uncertainty of a state 
localized in a quantum well, ~ h, to be much smaller 
than the total phase-space volume of the well. The latter 
can be estimated from the separatrix area, oc he/a ^ h, 
cf. second inset of Fig. [3]ja). Since in the semiclassical 
limit tunneling between the wells is exponentially sup- 
pressed, it is consistent to treat noise as local fluctuations 
in each well separately. 

Below the threshold, e < F, Eq. (El) imposes con- 



straints on the input field Bq. To be consistent with the 
linear description of fluctuations, we consider the quasi- 
linear limit of the classical response, a|Aop <C VF^ — e^. 
With this assumption, the maximum magnitude of the 
field in the cavity, Eq. (51 1 with Ob = — 7r/4, reads 



1^1 



2Fo 



(F-6)2 

while the amplified vacuum noise is 

£2/2 



F2 



(E4) 



(E5) 



The number of amplified vacuum photons inside the cav- 
ity is small at weak pumping but grows and passes the 



one-photon level at e = •\/2/3 F, and becomes dominant 
while approaching the parametric threshold. Using these 
estimates we extract from Eq. ( |E1[ ) the lower bound on 
the input signal. 



> 



2F2 ' 
e2 F-e 

4r2 F -I- e 



< 



2 e 1 

3<f <^ 



(E6) 



For very small pump strength, e ^ F, the constraint (E6 ) 
reduces to l-BoP ^ F, which is qualitatively similar to a 
high quality Duffing cavity, in which the resonant field 
fed by the input |i?oP ^ F achieves a large (classical) 
value, |Ao|2 ^ 1. The quasilinear approximation in this 
case, a <^ a\Ao\'^ <C F is valid as soon as a ^ F, and it 
imposes an upper bound on the input, |i?oP ^ F2/2a. 

Close to the threshold, the amplified signal grows with 
e more rapidly than the noise, and remains dominant 
at practically all input signals. This regime persists until 
the nonlinear effect breaks the quasilinear approximation 
at al^oP ~ VF2 — e2, and the signal amplitude satu- 
rates. The corresponding constraint on the input reads. 



\Bo\ 



F 



/2a 



5/2 



F-e<F. 



(E7) 



In terms of e, the upper bound for this regime is given 
by the condition 



(E8) 



For experimentally relevant cavity parameters, a/T < 
1/10, our estimates for the relative noise strength are 
therefore valid up to e ^ 0.95F. 
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